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ABSTRACT 

The S-stars, discovered in the vicinity of the massive black hole (MBH) in the Galactic center (GC), 
are anticipated to provide unique dynamical constraints on the MBH spin and metric, in addition to 
the mass. In this paper, we develop a fast full general relativistic method to simultaneously constrain 
the MBH mass, spin, and spin direction by considering both the orbital motion of a star close to 
the GG MBH and the propagation of photons from the star to a distant observer. Based on the 
current observations and dynamical model predictions, we assume six example stars with different 
semimajor axes (oorb) and eccentricities (eorb) and numerically calculate their projected trajectories 
in the sky plane and redshift curves. Two of those stars are set to have orbital configurations similar 
to that of S0-2/S2 and SO-102. We find that the spin-induced effects on the projected trajectory 
and redshift curve of a given star, including the leading term by the Lense-Thirring precession and 
the frame dragging, and the high-order precession due to the quadruple moment, depend on both 
the absolute value and the direction of the spin. The maximum values of the spin-induced position 
displacement and the redshift differences of the star over a full orbit may differ by a factor of several 
to more than one order of magnitude for two cases with significantly different spin directions. The 
dependence patterns of the position displacements and redshift differences on the spin direction are 
different, and thus the position and the redshift data are complementary for constraining the MBH 
spin and its direction. Adopting the Markov Chain Monte Carlo fitting technique, we illustrate that 
the spin of the GC MBH is likely to be well constrained by using the motion of S0-2/S2 over a period 
of ~ 45 years if the spin is close to one and if the astrometric and spectroscopic precisions can be 
as high as (cTp, az) ~ (10/ras, lkms“^), which is expected to be realized by future facilities like the 
GRAVITY on the Very Large Telescope Interferometer, the thirty meter telescope, and the European 
extremely large telescope. If Cp and az can be further improved by a factor of several, the MBH spin 
can be well constrained by monitoring S0-2/S2 over a period of ^ 15 years. In the mean time, the 
distance from the Sun to the GC and the MBH mass can also be constrained to an unprecedented 
accuracy (0.01%-0.1%). If there exists a star with a semimajor axis that is a few times smaller than, 
and eccentricity larger, than those of S0-2/S2, the MBH spin and its direction can be constrained 
with high accuracy over a period of < 10 year by future facilities, even if the spin is only moderately 
large. Our results suggest that long-term monitoring of the motions of stars in the vicinity of the GC 
MBH by the next generation facilities is likely to provide a dynamical test, for the first time, to the 
spin and metric of the GC MBH. 

Subject headings: Black hole-physics - gravitation - Galaxy: center ~ Galaxy: nucleus - relativistic 
processes - stars: kinematics and dynamics 


1. INTRODUCTION 

The Kerr black hole (|Kerdll9^ is one of the most sim¬ 
ple and elegant solutions to the Einstein field equation 
of general relativity (GR). It is widely accepted that all 
astrophysical black holes (BHs), if existing, can be de¬ 
scribed by the Kerr metric with only two parameters, 
i.e., the mass and the spin (no-hair theorem). Various 
lines of evidences for the existence of BHs have been 
accumulated in the past several decades, mainly based 
on dynam ical measurements of the masses of these ob¬ 
jects (e.g..lMagor rian et a nil998llKormendv fc Holl201,1I 
iNaravan fc McGlintockll201,'l[l . One of the strongest evi¬ 
dences is provided by the long-term monitoring of the 
motions of stars in the Galactic center (GC), which 
suggests, almost exclusively, the existence of a massive 
black hole (MBH; with mass ~ 4 x IO^Mq) in the GC 


(see Ghez et al.l[2008l iGillessen et al.l[20(i^ iMever et al.l 
12012 1. However, it is still not clear whether these objects 
can be fully described by the Kerr metric or not. 

A crucial step to check whether an astrophysical BH 
is described by the Kerr metric is to accurately con¬ 
strain/measure the BH spin through spin-induced GR 
effects. Recently, some progress has been made in mea¬ 
suring th e spins of both stellar-m ass BHs in the Milky 
Way li.e.. lMcClintock et al.l[2r)Tll and references therein) 
and MBHs in t he centers of active galactic nuclei (e.g., 
iRevnoldi I20I3L and references therein). These mea¬ 
surements are obtained by modeling the intrinsic X-ray 
continuum and its reflection components (e.g., the Fe 
Ka line and the reflection continuum) from each BH- 
accretion disk system, which may suffer from various un¬ 
certainties in accretion disk models. Therefore, indepen¬ 
dent determinations of the spins of BHs are of fundamen- 
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tal importance for testing the no-hair theorem, GR, and 
a deep understanding of the space and time. 

The MBH system in the GC, massive and in close prox¬ 
imity to us, provides a unique laboratory for testing the 
no-hair theorem and the GR due to the following reasons. 
First, the angular size of the expected shadow/image of 
the GG MBH is the largest one among those of MBHs, 
which makes the GG M BH the best target for the event 
horizon telescope fe.g.. iDoeleman et al.ll2009tl . Second, 


the GC MBH is the only MBH which has stars in its 
vicinity that can be detected individually. Long-term 
monitoring of the motions of those stars in the imme¬ 
diate vicinity of the MBH offers a novel opportunity to 
probe various GR effects, including the periastron ad¬ 
vancement, the Lense-Thirri n g prec e ssion, and the frame 


200(1: IR.ubilar h EckartI 120011: 

Weinberff et al.l 20051: 

Preto & Sahal 120091: lAnerelil et al 

20I0I: lAnffelil & Saha 

20IOI: IMerritt et al.l I20I0I: lAnsrelil & Sahal I2011b llorio 

20IIallbr). Gurrentlv, the dete 

cted S-stars rotating 
the smallest semima- 
32, which have semi- 
U and pericenter dis- 

around the GC MBH that have 
jor axes are S0-2/S2 and S0-I( 
major axes of 980 AU and 850 A 
tances of 120 AU and 270 AU, resp 

ectivelv (ISchddel et al. 

20I2:IGhez et al.ll2008HGillessen et al.ll2009HMever et al. 


201211 . It is also anticipated that there are some stars 
(fainter than the S-stars) and p ulsars hiding within the 
orbit of S0-2/S2 and SO-102 fsee lZhang et al.ll201,1l . 12013 . 
and references therein). The GR effects, e.g., the Lense- 
Thirring precession and the frame dragging, which are 
inherited from the MBH spin, shoul d hide in the ap- 
parent trajectories of these stars (e. g., 


and in the redshift curve of the star (lAngelil et al.ll201C 
however, they are almost unreachable by current facili¬ 
ties. 

A number of next-generation facilities are currently in 
progress, such as the GRAVITY on the Very Large Tele¬ 
scope Interferometer (VLTI), the thirty meter telescope 
(TMT), and the European extremely large telescope (E- 
ELT). With these facilities, not only will there be great 
advancement in the precisions of the astrometric and 
spectroscopic measurements, but also some stars, fainter 
and closer to the cent ral MBH than the known S-stars, 
may be revealed fe.g.. lZhang et al.ll2013[l . It is promis¬ 
ing to dynamically constrain the MBH spin and metric 
by using the relativistic orbital motions of those stars. 

Various spin-induced GR effects hav e been in¬ 
vestig a ted by many authors (e.g., I.Taroszvnskil 
1998t [Rubilar fc EckartI I200H lAngdlil et al.l 120101: 


Angdlil &: Sahal201^lMerritt et al.lf2010HAngelil &: Sahal 

201lHlorioll2nilallbfl , however, most of those studies are 
based on perturbative approximations rather than full 
general relativistic calculations and do not consider how 
tight the constraint on the MBH spin could be obtained. 
In this paper, we develop a full general relativistic 
method to investigate the spin-induced GR effect by 
considering both the relativistic orbital motion of a star 
in the vicinity of the GG MBH and the propagation of 
photons from the star to a distant observer. Using this 
method and the Markov Chain Monte Carlo (MCMC) 
fitting scheme, we investigate how a constraint on the 
MBH spin and how accurate the constraint can be 
obtained by using the orbital motion of an example 
star, e.g., S0-2/S2, for a given set of astrometric and 


spectroscopic accuracies of an instrument or a telescope. 

This paper is organized as follows. By adopting the 
Kerr metric to describe the curved spacetime around the 
MBH, we introduce the motion equations for a star ro¬ 
tating around the MBH and for photons propagating in 
the spacetime in Section [2j We introduce the numer¬ 
ical method to solve the GR motion equations for the 
star and the ray-tracing technique to trace back those 
photons propagating from the star to a distant observer 
in Sections [3] and SI In Section 5, we summarize the 
main results of our calculations on the orbital motions of 
some example stars and the propagation of photons from 
those stars to the distant observer. Two of those example 
stars have similar orbital configurations as those of SO- 
2/S2 and S0-I02; and the others are set according to the 
predictions on the probability distributions of the semi¬ 
major axis and eccentricity of a star, which is expected to 
be the closest one to the GC MBH as predicted by some 
dynamical models. In Section [51 we describe the pro¬ 
cedures of an MCMC fitting to the mock observations 
of the example stars, assuming astrometric and spectro¬ 
scopic accuracies of the future observations. By doing 
this, we further investigate how accurate the constraint 
on the MBH spin and metric can be obtained by using 
the relativistic orbital motions of a star, e.g., S0-2/S2, 
with the next generation facilities. Discussions on some 
related issues are given in Section [7] and conclusions are 
summarized in Section |8l 


2. KERR METRIC AND MOTION EQUATIONS 

The curved spacetime around an as trophysical BH may 
be described by the Kerr metric (iKer il 1 1 96311 . which 
can be expressed in the Boyer-Lindquist coordinates 
{t, r, 9, (j)) as 

ds2 ^ + (I) 

where 

' = SA/A, 

= Asin^e/E, 
e2^i = S/A, 

„2/i2 _ V 

> I = tr/A. (2) 

S = r2 -I- cos^ 9, 

A = r2 — 2r -I- a2, 

A = (r2 -I- — a2Asin20, 


and a is the dimensionles s spin parameter of the MBH 
(jBover fc Lindoui^ II967I 1. For simplicity, we set G = 
c = M, = Tg = GM^jt? = I above, G, c, M, and Vg are 
the gravitational constant, the speed of light, the MBH 
mass, and the gravitational radius, respectively. 

In the Kerr metric, the motion of a particle, a star, 
or a photon, is controlled by the following equations 
(iBardeen et al.l[T972HMisner et al.llI973HChandrasekhad 

fl^ 


Sr = ±\/R, (3) 

T.9 = ±y/Q, (4) 

E(/ = —a-I-A/sin^ 0-I-oT/A, (5) 

St = —sin2 0-I-oA-I-(r^-|-a^)r/A, (6) 


where overdot represents an ordinary derivative with re¬ 
spect to an affine parameter r, and T = r^-|-a^ —Aa. The 
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quantities R and 0 in Equations ([3]) and O are given 

by 

R={1- ^2)^4 2^2^3 _ ^2) _ ^2 _ ^2]^2 

+2[(a-Af+ g>-aV, (7) 

and 

0 = g 2 _ [^2^^2 _ ^2^ ^^^^2 ^^^2 

respectively, where ^ = m/E = 0 for photons), A = 
LzjE, and = Q/E"^. The constants of motion are 
the rest mass m, the energy at infinity E, the azimuthal 
angular momentum Lz, and the Carter’s constant Q. 

For a star in the immediate vicinity of an MBH that 
is monitored by a distant observer, the GR effects en¬ 
coded in the measurements can be divided into two parts. 
First, the orbital motion of the star rotating around the 
MBH is affected by the GR effects, such as the peri- 
astron advancement, the Lense-Thirring precession, and 
the frame dragging. Second, the propagation of photons 
from the star to the distant observer is also affected by 
the GR effects, including the gravitational redshift for 
those photons to climb out of the MBH potential, and 
the gravitational bending of the photon trajectories if the 
distance of those photons to the central MBH is small. 
Therefore, it is necessary to consider both the relativistic 
motion of the star and the propagation of photons from 
the star to the distance observer, simultaneously. Below, 
we introduce the full general relativistic numerical meth¬ 
ods to calculate the relativistic orbits of stars and the 
trajectories of photons in the Kerr metric, respectively. 

3. RELATIVISTIC ORBITS OF STARS 

The orbital motion of a star in the vicinity of the 
central MBH, with a fixed mass (M,), spin (a) and 
its directionQ is determined by equations ©-([S]) once 
the initial conditions are known. The initial condi¬ 
tions of the star can be given by its position and tetrad 
velocity at a time t*.o in the Boyer-Lindquist coor¬ 
dinates, i.e., r*^o = (G.o, 'r*.o, ^**. 0 , and u*_o = 
(it* Q, it*Q, It* oi ^*o)' Hereafter, t*_o is set to be zero 
for simplicity, if not otherwise stated. At any given time 
t*, the three-velocity v* = (i!*,u*,i'f) in the space of 
the local non-rotating rest frame (LNRF) associated with 
(r*,0*,(()*) can be transformed from the tetrad velocity 
u* = as (s ee iBardeen et al.l 

119721 : Il3ajdeenlll973l : iMisner et al.lll973l) 


vl = 

(9) 


(10) 

= {ut/ul - uj) 

(11) 


The tetrad velocity u* can be reversely transformed from 
the three-velocity v* in the LNRF frame. 

If the position and the three-velocity of a star in the 
LNRF frame are provided, the values of ^ and A can be 
derived as 

^ = l/[j{e''+coe^vt)] ( 12 ) 

^ The MBH mass, spin, and its direction are assumed to be fixed 
values because the time duration of the orbital motion considered 
in this study is short and no more than several decades, and thus 
there should be no significant evolution in the MBH mass and spin 
vector. 


and 

X = Avt/i2arvt + r), (13) 

respectively, where 

7=(1-v*-v*)-H2, (14) 

and 

F = [A( A - sin^ 6) + sin^ 9] j E^A sin 9 . (15) 

Substituting ^ and A into Equations © and m , we 
have 

= E (yuf)^ -f cos^ 9 [a2(^2 _ i) AVsin^ 9] . (16) 

We define two pseudo-Cartesian coordinate systems for 
the convenience of relating the motion of a star in the 
vicinity of the MBH with the Boyer-Lindquist coordi¬ 
nates (r, 9, (j)) to that seen by a distant observer. The 
first pseudo-Cartesian coordinate system x'y'z' is defined 
in the rest frame of the distant observer (denoted as the 
observer’s frame hereafter). The direction of the z'-axis 
is pointing from the MBH to the distant observer, and 
the x'y' plane represents the sky plane of the observer. 
The second pseudo-Cartesian coordinate system xyz is 
defined to relate the Boyer-Lindquist coordinate (r, 9, (f>) 
to an orthogonal coordinate system (x, y, z) with z rep¬ 
resenting th e direction of the s pin vector (see a similar 
definition in lAngelil et al.l[2010ll . Because the star is usu¬ 
ally far away from the event horizon of the central MBH, 
the coordinates of the star in the xyz frame may be ap¬ 
proximated by 

X = rsin0cos()), 

y = rsin^sinc^, (17) 

z = rcos9. 

In this xyz frame, the direction of the z-axis represents 
the spin direction of the MBH, the y-axis represents the 
intersection line of the equatorial plane of the MBH with 
the sky plane of the distant observer. The LNRF frame 
is rotating with an angular velocity of uj relative to the 
distant observer’s frame. However, uj is relatively small 
since the star is usually far away from the event hori¬ 
zon of the MBH. Therefore, the rotation of the LNRF 
frame with respect to the distant observer’s frame can 
be neglected when relating the initial conditions to the 
Newtonian approximations. 

Figure [T] shows the configuration of the stellar-MBH 
system and the two pseudo-Cartesian coordinate systems 
defined above. As seen from Figure [U the spin direction 
z is determined by two angles in the observer’s frame 
x'y'z', i.e., the angle (i) between zand z', and the angle 
(e) between the projection of .?on the x'y' plane and —x'. 
The angle between the y axis and the y' axis is the same 
as e. Therefore, a vector defined in the observer’s frame 
can be approximately transformed to the LNRF frame 
(or the MBH’s frame) by first rotating the x'y'z' frame 
around the z' axis counter-clockwise by an angle of e and 
then rotating it around the new y' axis clockwise by an 
angle of i. For example, a vector n' in the observer’s 
frame can be transformed to that in the LNRF frame n 
by n = Mn', where M is a rotation matrix given by 

( cos* cose cos i sine sin A 

— sine cose 0 J . (18) 

— sin i cos e — sin i sin e cos i/ 
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To the observer 



Fig. 1. — Schematic diagram for the star-MBH system and the coordinate systems. Right panel; a pseudo-Cartesian coordinate {x',y', z') 
is defined in the rest frame of a distant observer (i.e., the observer’s frame), where the x'y' plane represents the sky plane of the observer 
and is taken as the reference plane, the z' axis represents the line of sight, and the x' axis is taken as the reference direction on the sky plane. 
Another pseudo-Cartesian coordinate xyz is defined to relate the Boyer-Lindquist coordinates to orthogonal coordinates (x,y,z), 

where the z axis represents the spin direction of the MBH and the y axis represents the intersection line of the MBH equatorial plane with 
the observer’s sky plane. The direction of the MBH spin in the observer’s frame is therefore described by two angles, i.e., the angle the z 
axis and the z' axis (f) and the angle between the projection vector of the x axis on the x'y' plane and the x' axis (e). The directions of 
the R.A. {—y') and Dec. {—x') are marked in this figure. Left panel: the orbit of the star may be approximated as a Newtonian ellipse in 
the observer’s frame, especially when the semimajor axis of the star is large, which is described by six orbital elements, i.e., the semimajor 
axis Uorbi eccentricity e, the longitude of ascending node Q', argument of periapsis T', true anomaly v' and inclination I', respectively. 


In the Newtonian case, the orbital motion of a star 
rotating around a massive object is determined by six 
orbital elements, i.e., the semimajor axis Oorb, the eccen¬ 
tricity Corb, the longitude of the ascending node fl', the 
argument of the pericenter T', the true anomaly v\ and 
the orbital inclination angle I' (see Figure [1]). For con¬ 
venience in comparison with the Newtonian orbits cur¬ 
rently determined for a few GC S-stars with the smallest 
semimajor axis, in our following simulations of the or¬ 
bital motions of some (example) stars, we set their initial 
conditions by fixing the six orbital elements at a given 
moment in the distant observer’s frameO Then, the val¬ 
ues for the four-position r,t,o f^d the tetrad-velocity 
of a star can be approximately obtained by the following 
procedures. 

1. The initial position {x'^ o,yi,o, z'^ q) and the three 

velocity of the star in the distant 

observer’s frame are first obtained from the six or¬ 
bital elements initially set. 

2. Transforming the position and the velocity of the 

star in the distant observer’s frame into those in 
the LNRF frame, i.e., ry,o = (a;*,o, y*,o, •z*,o) and 
fti.o = by the rotation given in 

Equation la, and then transforming these two 
vectors from that using the xyz coordinates into 
that using the r9(j) coordinates. 

3. Transforming r*_o and to the four-position r*-_o 
and the tetrad-velocity u* g in the Boyer-Lindquist 
coordinates according to Equations dlD-dlll)- 

^ The orbital elements set here are used in a simple way to gen¬ 
erate the initial conditions of a star moving in the Boyer-Lindquist 
coordinates. The elliptic orbit defined by these orbital elements 
is only taken as a Newtonian approximation to the real orbit of 
a star in the Boyer-Lindquist coordinates. These orbital elements 
may lose their original meaning in the curved spacetime around a 
Kerr MBH. 


Once the initial conditions are given for a star, the mo¬ 
tion parameters A, and q can be obtained by Equations 
(O, dSl), and m, respectively. With those constants 
of motion, the orbital motion of a star, in principle, can 
be obtained by integrating the motion Equations (I3|)-(l6]). 
However, it is required to frequently judge the sign “±” 
at the left side of the motion Equations ([3]) and ([4]) in the 
numerical integrations. To avoid this complication, one 
may use the following equations to replace the motion 
Equations (HI) and ([4]). 


<1 

II 

w 

(19) 

^0=Pe, 

(20) 

1 dR 50 

“ 2A 5r “ ^ 50 ’ 

(21) 

150 

2 50■ 

(22) 


We use the code DORPI5 based o n the explicit fifth 
(four t h)-order Ruirge Ku tta method (|Dormand fc: Princ^ 
119801 : iHairer et al.l ll993D to integrate the motion Equa¬ 
tions (1I1|-(E21) and dS]) and ([HI) to obtain the orbit motion 
of a star. We set the relative integration errors to be 
< 10“^^ for all the position and momentum quantities 
(r, 0, (j), pr, pg, p^) in those equations, which are suffi¬ 
cient for the convergence of the numerical results and 
the required accuracy in this study. 

3.1. Orbits of Example Stars 

Observations in the past two decades have revealed a 
number of GC S-stars rotating around the MBH with 
semimajor axes in the range of ^ 4000-800 AU. Among 
them, the two wi th the smallest semimajor axes are SO- 


2/S2 and SO-102 

Ghez et alJl2008l: IGillessen et alJl2009l: 

IMever et al.l 12012 

). These stars have been very useful 


in determining the MBH mass, and it is expected that 
S0-2/S2 (or a star within the orbit of S0-2/S2) will be 
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Fig. 2.— Orbital distribution of some detected GC S-stars and the possibly existing star that is closest to the GC MBH in the semimajor 
axis versus eccentricity plane (left panel) and in the semimajor axis versus the distance to the pericenter plane (right panel). The tick labels 
at the right side of each panel and the top side of the right panel are in units of gravitational radius rg = GM,/c^, while the tick labels at 
the left side of each panel and the bottom side of the right panel are measured with the astronomical unit AU. The solid circles represent 
those currently detected GC S-stars within a distance of 400 0 AU from the GC MBH. S0-2/S2 and SO-102 are the two with the smallest 
semimajor axes among the GC S-stars fsee lMever et akll2012| j. The probabilities of a predicted star with mass in the range 1 — 7Mq that 
is closest to the MBH within t he yellow and the m agenta solid (dotted) curves are 85% (98%) and 66% (98%), resulting from the model 
“Disk-IM2” and “Disk-IMO” in IZhang et al.l II2013I1 . respectively, by assuming that the GC S-stars are the captured components of stellar 
binaries injected into the vicinity of the MBH. The cyan dashed lines give the gravitational radiation timescale of a star with mass SMq 
rotating around the MBH: 10®, fO^, and fO® years from top to bottom in the left panel (or from right to left in the right panel). The 
pericenter distances of the stars, if existing, in the shaded region are smaller than the tidal disruption radii tor a star with mass 3Mq, and 
thus stars almost cannot exist in this region. 


an important dynamical probe to the GR effects of the 
GC MBH. It is also anticipated that there are some stars 
fainter than the GC S-stars existing within the orbits of 
S0-2/S2 and SO-102, which could b e even better probes 
for the GR effects. iZhang et al.l (I2013D predicted the 
probability density distribution in the semimajor axis 
versus eccentricity plane of such stars (with mass in the 
range of 1 — 7Mq), by assuming that these stars and the 
GC S-stars are the captured components of stellar bina¬ 
ries that are tidally broken up in the vicinity of the GC 
MBH. 

Figure [2] shows both the orbital distribution of the 
detected GC S-stars (solid circles) and the probabil¬ 
ity density distribution of the expected innermost stai 0 
(color contours) in the semimajor axis versus eccentric¬ 
ity plane (panel (a)) and the semimajor axis versus peri¬ 
center distance plane (panel (b)). The magenta (yel¬ 
low) solid and dotted contours represent a probability 
for the expected innermost star with mass in the range 
of 1 — 7Mq of 85% and 98% (or 66 % and 98%), re- 
sulted from t he “D isk-IMO” (or “Disk-IM2”) model in 
IZhang et al.l ()201.8[) . respectively. The semimajor axis 
and eccentricity of the expected innermost star are in 
the range of ~ 100 — 1000 AU and 0 — 0.99, respec¬ 
tively. The dotted lines in Figure [2] give the gravita¬ 
tional radiation timescale of a star with mass of 3Mq 
rotating around the MBH with mass 4 x I0 ®Mq: 10®, 
10^' 10® yea rs from top to bott om (Equation (5.6) i n 
iPetersI [T964I or Equation (39) in I Yu fc Tremaingl2003D . 
A star should be tidally disrupted if it approaches the 

MBH within a distance of rtid = r* ~ 

^ Here the “innermost” star means the star with the smallest 
semimajor axis. 


1 /3 

(ri'^M,/m^) (TO*/M 0 )°'^^i? 0 , where to* is the mass 

of the star and rj is set t o be 2.21 for a homogeneous , 
incompressible body Isee iMagorrian fc Tremaingll999D . 
The upper (right) boundary of the shaded region in the 
left (right) panel of Figure [2] represents the tidal radius 
of stars (rtid) with mass 3 M 0 rotating around the MBH. 
In principle, stars with masses > 3 M 0 cannot exist in 
the shaded region. The probability of existing stars in 
the region between the top (right) dashed line and the 
upper (or right) boundary of the shaded region in the 
left (or right) panel is substantially suppressed because 
of the rapid orbital decay of a star in this region due to 
gravitational wave radiation. 

According to the current observations and theoretical 
expectations, we select six example stars as listed in Ta¬ 
ble [1] Two of the example stars have similar orbital 
properties as those of S0-2/S2 and SO-102, respectively, 
which have the smalles t seinimajor axes among the de¬ 
tected GG S-stars (see Ghez et al.ll2008l : iGillessen et al.l 
120091 : [Mever et aI1l2012 ). The other four example stars, 
i.e., Ea, Eb, Ec, and Ed (see Table [T|), are set according 
to the probability distribution of the expected innermost 
star shown in Figure [2l The semimajor axes, eccentric¬ 
ities, and pericenter distances of the example stars Ea, 
Eb, Ec, and Ed are (300 AU, 0.88, 36 AU), (300 AU, 0.98, 
6 AU), (80AU, 0.88, 9.6AU), and (80AU, 0.3, 56 AU), 
respectively. The probability for the existence of a star 
like the example star Ea is probably high. There are 
some chances for the existence of a star like the example 
stars Eb, Ec, or Ed in the vicinity of the MBH, though 
the chances are not very high. Note that the semimajor 
axis and eccentricity listed there are only taken as the 
Newtonian approximation to the orbit of each example 
star, which are useful for setting the initial conditions. 
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TABLE 1 

Orbital parameters for example stars 


Name 


*^orb 


■ Sorb 

I 

n' 

'Y'/ 

'-pd 

-^0 

a 



AU“ 

r ^ 

mas 

S0-2/S2 

984 

24949 

123 

0.88 

135° 

225° 

63° 

2.32 

0.99 

45° 

200 ° 

SO-102 

848 

21500 

106 

0.68 

151° 

175° 

185° 

9.5 

0.99 

45° 

180° 

Ea 

300 

7606 

37.5 

0.88 

45° 

0 ° 

0 ° 

0 

0.99 

45° 

180° 

Eb 

300 

7606 

37.5 

0.98 

45° 

0 ° 

0 ° 

0 

0.99 

45° 

180° 

Ec 

80 

2028 

10 

0.88 

45° 

0 ° 

0 ° 

0 

0.99 

45° 

180° 

Ed 

80 

2028 

10 

0.30 

45° 

0 ° 

0 ° 

0 

0.99 

45° 

180° 


Note. — “In units of AU. ** In units of the gravitational radius Vg = and the MBH mass M, is assumed to be 4 X 10® Mq. “ 

In units of mas by assuming the distance from the MBH to sun i?GC = 8 kpc. Time of pericentric passage, with respect to the year 2000 

in order to obtain mock observations for S0-2/S2, SO-102, and other possibly existing stars by future facilities. “ The MBH spin direction 
is set for the case for S0-2/S2 in order to have an intermediate spin-induced effect; and it is set for the cases of other stars arbitrary to the 
values listed in the table. 


4. PHOTON PROPAGATION: RAY-TRACING 

We consider a photon propagating from a star in the 
vicinity of an MBH to a distant observer in this section. 
The position of the star at each moment (t*, r*, 0*, 0*) is 
set by its orbital motion as obtained above in Section [3] 
once the initial conditions are fixed. The position of the 
distant observer is fixed at (to, 0o, f/'o) = (i?GC,*,0°), 
where i?GC is the distance from the Sun to the GC 
and set to be 8 kpc if not otherwise stated. The pho¬ 
ton trajectory may be bended due to the curved space- 
time of the MBH, and the energy of the photon when it 
is received by the distant observer is different from its 
original value emitting from the star due to the grav¬ 
itational redshift and the Doppler shift. To account 
for those effects, we use the ray-tracing technique to 
trace back photons from the distant observer to the 
star in the Boyer-Lindquist coordinate system. The 
ray-tracing technique is detaile d in Appendix 1X1 where 
both Jacobian elliptic functio ns (|Bvrd fc FriedmaJll954 
lAbramowitz fc Stegunlll97^ and the Gauss-Kronrod in¬ 

tegration scheme le.g.. see lRauch fc Blandfordlll994l l are 
adopted to integrate the motion equations. 

The apparent position of the star in the distant ob¬ 
server’s sky plane at a given time to is described by two 
impact parameters, (a, /3), which map the position of 
the star (t*, r*, 0*, </>*) and are measured relative to the 
center of the MBH. The impact parameter a is the ap¬ 
parent displacement of the image perpendicular to the 
projected axis of symmetry of the MBH, and /3 is the 
apparent displacement parallel to the projected axis of 
symmetry in the sense of the angular momentum of the 
MBH. For photons, ^ = 0, the values of the other two 
constants of motion, A and in Equations o and (|8|), 
can be determined by a and /? as 

A = —asinf, (23) 


and check whether the photon can reach the surface of 
the star at a given moment. 

According to the motion Equations dS]) and , we 
have 



(25) 


For a photon propagating from the star to the distant ob¬ 
server, the 0-integral at the right side of the above equa¬ 
tion can be obtained by i ntegr ating over 0 from 0o to 0* 
(for details, see Appendix I A. II) . When 0 = 0*, the radial 
position of the photon is denoted as rhit- Therefore, the 
r-integral at the left side can be obta ined by integrating 
over r from To to rhu (see Appendix lA.21) . By equating 
the r-integral with the 0-int egral , we can then get the so¬ 
lution of Thit (see Appendix IA.3I) . Once rhit is obtained, 
one can subsequently obtain the azimuthal position of 
the photon, c/ihit, when 0 = 0* and r = rhit, and the time 
needed dthit = thit — to for the photon propagati ng fr om 
(thit,’’hit,(Ahit) to the observer (see Appendix IA.4p . 

The trajectory of a photon with arbitrary (a, /3) from 
the distant observer may not hit the surface of the star 
at any time t* when 0 = 0*. We adopt the following pro¬ 
cedure to judge whether a photon can reach the surface 
of the star. In the pseudo-Cartesian coordinate system 
xyz, we have i?hit = (iEhit, 0hit, ^^hit) and 


a;hit = ’’hit sin 0* cos ^hit, 

i/hit = ’’hit sin 0* sin <(>hit, (26) 

-^hit — ’’hit COS 0*, 


according to Equation ini). The position of the star at 
any given time t* can be described as i?* = (a;*,y*,z*), 
which are obtained similarly as i?hit- The distance of the 
photon at 0 = 0* to the star in the xyz system is defined 
as 


= ^32 _i_ ^^2 _ ^2^ ^pg2 ^ 24 ) 


d = 


.^hit 


For details, see Cunningham k, BardeenI (1197311 . 
iChandrasekhad (|1983 1. and Karas et al.l ( 199^ . The 
star may have more than one image because of the 
gravitational leasing effect. For the majority of the cases 
investigated in this study, however, only one image is 
dominant because the star is normally far away from the 
Einstein radius of the central MBH and the leasing effect 
is not significant. We may trace a photon with given 
impact parameters (a, (3) from the distant observer to 
the vicinity of the MBH by solving Equations (Ill)-®, 


= V(aihit - a;*)2 -|- (j/^it - y*Y + (^hit - -2*)^- (27) 

We assume that a photon trajectory hits on the surface 
of the star \i d < and ^ = 10“^ — 10“®13 In this case, 

^ For a star with mass 3 Mq, its radius is ~ 1.7x 10® 0.3rg 

if M» = 4 X 10®M 0 . The setting of ^ = 10“® — 10“® enables suf¬ 
ficient accuracy for relating the position of an example star to its 
image in the observer’s sky plane, as we have d < 0.05rg, substan¬ 
tially less than the star radius, even for the example star with the 
largest semimajor axis. 
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(a, P) represent the apparent position of the image of the 
star at (t*, r*, 0*, (/>*) in the observer’s sky plane, which 
are re-denoted as (a*,/3*). The image (a*,/?*) should 
be detected at a time to = Therefore, the 

orbital motion of stars in the vicinity of the MBH rep¬ 
resent by the change of a* and /3* in a time sequence 
of to- In the meantime, the four-momentum of a pho¬ 
ton with motion parameters A and q at any position, 
^•§■5 Phit — (Pt,hit5 Pr,hit;P0,hit5 P0,hit) (r'hit; ; ^hit) 

and Po = (pt.o,Pr.o,Pe,o,P0,o) at (i?GC,*,0°), can be ob¬ 
tained by 

Pt = -E, 


Pr = eeVH/ a, 

pg = ±eVq, 

, = \E, 


(28) 


in the Boyer-Lindquist coordinates. We can then obtain 
the redshift of the photon as 

^ ^ ^ ^ _Phn^ _ ^ ^ 29 ) 

Po ■ Uo Eo 

Here Uq = (—1,0,0, 0) and u* are the four-velocities of 
the distant observer and the star, respectively, Eo is the 
energy of the photon received by the observer. According 
to Equation the redshift curve can be obtained for 
each star. This redshift can be measured by the shift of 
the center(s) of emission or absorption line(s). 

In order to perform the ray-tracing calculations effi¬ 
ciently, the integrations are started from the point at 
(10®rg, i, 0) rather than (i?GC) h 0)- This setting does not 
affect the estimations of the impact parameters a and (3 
since the bending of light is not important anymore at 
a distance larger than lO^r-g from the MBH. The final 
(mock) observations of the apparent positions of a star 
on the observer’s sky plane should be given by ar^/Rcc 
and /3rg/i?GC- 


5. APPARENT ORBITS AND REDSHIFT CURVES OF 
EXAMPLE STARS 

The apparent position of a star in the observer’s sky 
and the redshift of the star at any given time can be ob¬ 
tained according to the procedure described in Sections [3] 
and m We illustrate the effects of the MBH spin and its 
direction on the apparent orbits and redshift curves by 
using those example stars listed in Table[T]in this Section. 


5.1. Spin-induced Apparent Position Displacements 

Figure 13] shows the trajectory of each example star on 
the observer’s sky, which is described by the right as¬ 
cension (R.A.) and declination (Dec.) of the star with 
respect to the GC MBH (or Sgr A*). The initial set¬ 
tings for the orbital parameters of each star are listed in 
Table |T] As seen from Figure |31 the orbit of each star 
precesses mainly due to the advancement of periapsis; Eb 
has the most significant advance precession, while SO-102 
has the least advance precession. The advance precession 
of periapsis of a star per orbit is 


m=37r 


2r„ 


aorb(l ^orb) 


-0.19° 


4 X 10® Mq 

M, 


/10® AUA 1-0.882 


“"orb 


where Oorb and Corb are the semimajor axis and eccen- 
tricity of the star i n the Newtonian approximation (e.g., 
iMisner et aLlHOT^ . For those example stars listed in 
Table (H JT'g are 0.19°, 0.095°, 0.64°, 3.64°, 2.39°, and 
0.59°, respectively. Considering the projection of the 
orbital plane of each star to the observer’s sky plane, 
the apparent shifts of the apoapsides of the example 
stars after three full orbits are 3aorb(l + eorb)i5Tg(l — 
cos2T'sin2/')i/2 ~ 2.23, 0.77, 1.67, 10.0, 1.67, and 
0.29 mas for S0-2/S2, SO-102, Fa, Eb, Ec, and Ed, respec¬ 
tively (see Fig. 131). According to Figure |31 the advance¬ 
ment of periapsis can be accurately measured within a 
few orbits for all of the example stars, including SO- 
2/S2 and SO-102, if the astrometric precision can reach 
100 — lO^as scale. The spin-induced effects, including 
the Lense-Thirring precession and the frame dragging, 
cannot be clearly seen in Figure |3] because they are more 
than one order of magnitude smaller than the periap¬ 
sis advancement, and these effects are addressed in Fig¬ 
ures [31131 

If the MBH is non-spinning (a = 0), the trajectories 
of those example stars on the observer’s sky plane at a 
given moment may slightly deviate from those shown in 
Figure 131 for the case with a = 0.99. Figured! shows the 
evolution of this difference for each star over three full 
orbits. The difference in the apparent position is mainly 
due to the Lense-Thirring precession and the frame drag¬ 
ging. Figure dl shows the distance {S) between the appar¬ 
ent position of each star on the observer’s sky plane at 
each given moment for the case with a = 0.99 and that 
with a = 0. After the motion of one full orbit of each 
star, the differences at apoapsis are — 8.7, 1.6, 4.7, 92.5, 
11.1, and 2.6^as for S0-2/S2, SO-102, Fa, Eb, Ec, and 
Ed, respectively (see Fig. El). These values are roughly 
consistent with simple analytical estimates as illustrated 
below. This difference mounts up with the increasing 
number of orbits as shown in Figures Eland El 

The Lense-Thirring precession and the frame dragging 
cause the precession of the orbital plane of a star, which 
leads to a difference between the position of the star in 
the sky plane at apoapsis/periapsis for the case with a 
rapidly spinning BH and that with a non-spinning BH af¬ 
ter one full orbit. For perturbation theory, this difference 
is approximately given by 

<^apo/peri - aorb(l ± Corb) [(5^D'(1 - sin^ T' sin^ /')-h 

i 52 t'(i _ eos^ T' sin^ I') sin T'^ sin^ I'S^I' -\- 

2 cos I'Sn'ST' - 2 sin T' cos T' sin /' cos I'ST'SI' 

- 2sinT'cosT'sin/'5D'(5/']^/^ (31) 

Here, the signs “-I-” and are for the difference at 
apoapsis (5apo) and periapsis ((jperi), respectively, JD', 
5T', and 51' are the changes of D', T', and I' due 
to the Lense-Thirring precession over one full orbit, re¬ 
spectively. For the detailed derivation of Equations EB 
see Appendix [BJ According to Equation (ISTI) . we obtain 
<5apo = 8.6, 1.5, 4.4, 62.5, 8.4, and 0.7/ras; Jperi = 0.55, 
0.30, 0.28, 0.64, 0.55, and 0.40 ^as for the example stars 
S0-2/S2, SO-102, Fa, Eb, Ec, and Ed, respectively, if 
a = 0.99 for the rapidly spinning BH case. The 5apo 
values are roughly consistent with numerical results ob¬ 
tained above for most example stars, except that for Ed. 

We note here that the 5 at periapsides shown in Fig- 
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Fig. 3.— Evolution of the apparent position of each example star on the observer’s sky plane. The apparent position of an example 
star is described by the R.A. and Dec. Panels from left to right, top to bottom, represent the example stars S0-2/S2, SO-102, Ea, Eb, Ec, 
and Ed, respectively. The initial settings of the orbital parameters for these stars are listed in Table ^ Three full orbits are shown for 
each example star. Open and solid circles mark the locations of the apoapsides and periapsides, respectively. In each panel, the dashed 
and dotted lines represent the eccentric vectors obtained by approximating the star orbit to a Newtonian ellipse for the first and the third 
orbit, respectively. As seen from the figure, the apparent orbital precession is most significant 3.1 mas per orbit) for the example star 
Eb, which has the smallest pericenter distance. 


ure[5]is obtained in a way slightly different from the 5apo 
estimated from Equation m- The values of S are cal¬ 
culated for the projected distances at any given time be¬ 
tween the star rotating around a rapidly spinning MBH 
(a = 0.99) and a star with the same initial orbital ele¬ 
ments but rotating around a non-spinning MBH (a = 0); 
while the value of 5apo is estimated by considering the 
difference between the position differences of two adja¬ 
cent apoapsides of a star rotating around an MBH with 
a = 0.99 and those of a star with the same initial orbital 
elements rotating around a non-spinning MBH. As the 
orbital periods of these two stars are slightly different, 
the shift of S at two adjacent apoapsides shown in Fig- 
ure[5]should therefore be roughly the same as, but differ¬ 
ent from, the analytical estimates ^apo, and the difference 
is most significant for those with small eccentricity. 

For a consistency check between the numerical results 
obtained in this study and the analytical ones, as an ex¬ 
ample, we show the apoapsis shifts numerically estimated 
in two different ways for those stars with the same semi¬ 
major axis of 300 AU (or 80 AU) but various eccentricities 
after one full orbit motion in Figure |6l The other initial 
orbital elements of those stars are the same as those of 
Ea (or Ec). The top panel of Figure IH shows the numer¬ 
ical results obtained in the same way as ^apo, which are 
well consistent with the analytical estimates at e < 0.8, 
but higher than the analytical ones at e > 0.9. There 
are three reasons for the difference at high eccentrici¬ 


ties: (1) the Newtonian approximation to the orbit of a 
star for the initial conditions becomes inaccurate if the 
star is close to the MBH (e.g., for those orbits with ex¬ 
tremely high eccentricities); (2) for an Ea-like star with 
Corb higher than 0.98, its orbital precession can be larger 
than ~ 20 per orbit, therefore, the first order approx¬ 
imation (for the projection) to obtain Equation (IHTl) is 
inaccurate; (3) the contributions from the higher order 
precessions are ignored in the analytical estimates. As 
showed in EigurelU the combinations of the above three 
factors can almost lead to 20% to a factor of two dif¬ 
ference when Corb > 0.98. The bottom panel of Figured 
shows the numerical results of S obtained in the same way 
as those shown in Figure[5l which apparently differ signif¬ 
icantly from the analytical estimates at both Corb ^0-7 
and Borb > 0.9 because S at apoapsis and Japo are defined 
and obtained in a slightly different way. 

According to Figure EJ we conclude that, if the ac¬ 
curacy in determining the apparent positions of any of 
those example stars (except SO-102) on the sky plane 
can reach < 10 /xas, the spin of the GC MBH can be 
then constrained by htting the evolution of its positions 
over several orbits as demonstrated by a Bayesian fitting 
method in Section |6] (see also Tab. [2]). 

For a star with given initial orbital elements, the dis¬ 
placement of apoapsis after one full orbit due to the spin 
effects is proportional to the absolute value of the spin, 
and it also depends on the spin direction (see eas. lB5llB8| . 
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Fig. 4. — Difference between the apparent position of each example star on the observer’s sky plane at each given moment for the case 
with a rapidly spinning central MBH (a = 0.99) and that with a non-spinning MBH (a = 0). The differences in R.A. and Dec. are denoted 
as S R.A., and (5Dec., respectively. Panels from left to right, top to bottom, represent the example stars S0-2/S2, SO-102, Ea, Eb, Ec, and 
Ed, respectively. Three total orbits for each star are shown here and the first orbit of each star starts at (^R.A., 5Dec.)= (0, 0). Solid and 
open circles mark the periapsis and apoapsis passage points, respectively. This figure shows that the difference in the apparent position of 
each example star increases with increasing time because of the Lensing-Thirring precession and the frame dragging and other higher order 
GR effects induced by the MBH spin in the rapidly spinning MBH case. 


To illustrate the dependence of the spin-induced posi¬ 
tion displacement on the spin direction, the top panel 
of Figure [7] shows the evolution of <5 of Ea for three 
different spin directions. The MBH spin is fixed at 
a = 0.99. As seen from this figure, S is the most signif¬ 
icant when (i, e) = (72°, 90°), and the least significant 
when (i, e) = (159°, 63°). The difference in the spin di¬ 
rection may lead to more than an order of magnitude 
difference in the maximum position displacement over a 
full orbit, which suggests that the MBH spin may be eas¬ 
ier to be constrained by the relativistic motion of a given 
star if the spin vector is close to some special direction. 
We note here that Japo of Ea is the largest (^apo — 20//as 
for a = 0.99) if (i,e) = (72°, 90°) or (108°, 270°), and 
it is almost the smallest (Japo — 0.4/ias for a = 0.99) 
if (i,e) = (159°, 63°) or (21°, 243°). (These values are 
obtained from Equation [STJ ) 

The spin-induced (apoapsis) position displacement (<5) 
also depends on the orientation of the orbital plane of a 
star relative to the sky pla ne, if th e MBH spin and its 
direction are fixed (see eqs. IB5IIB8I1 . The bottom panel 
of Figure [7] illustrates this dependence by plotting the 
evolution of <5 for Ea and two other stars with the same 
seinimajor axis and eccentricity as those of Ea but with 
different orbital orientations. The MBH spin is fixed at 
a = 0.99 and the spin direction is (i, e) = (45°, 180°). As 
seen from the figure, the amplitude of S may be different 
by more than an order of magnitude and the evolution 


pattern of S can also be different if the orientation of 
the star orbital plane is different, simply because of the 
projection effect. 

The apparent positions of a star, rotating around the 
MBH, on the observer’s sky plane, may be significantly 
affected by the bending of light rays propagating from 
the star to the observer due to the relativistic potential 
of the MBH. Figure [8] shows the difference in the appar¬ 
ent position of each example star at each moment for the 
case with, and without, the consideration of the bending 
of light (obtained from the ray-tracing method detailed in 
AnnendixIAT). As seen from this figure, the position differ¬ 
ence can be as large as 10—20/ias, roughly on the order of 
the Einstein radius of the MI3H (^ 4GM,/c^ ~ 20yLtas), 
which must be accurately considered when constraining 
the MBH spin, as the spin-induced position displace¬ 
ments are mostly on the same order (see Fig [S]). The 
position displacements caused by the bending of light 
are almost the same for the case with a rapidly spinning 
MBH and the case with a non-spinning MBH, and thus 
cannot be distinguished from each other (see the solid 
lines for cases with a = 0.99 in Figure [Si the results 
obtained for cases with a = 0 coincide with those solid 
lines, and the differences between these two cases are on 
the order of 0.1/ras). These differences, though small, are 
roughly on the same order as the effects induced by the 
quadruple moment of the MBH for S0-2/S2. Therefore, 
it is important to include them for future high-precision 
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Fig. 5. — Distance between the apparent position of each example star on the observer’s sky plane at each given moment tor the case with a 
rapidly spinning central MBH (a = 0.99) and that with a non-spinning MBH (a = 0). The distance is defined as <5 = ^(5R.A.)2 + (5Dec)^, 
and <5R.A. and (SDec. are shown in Figure|4] Panels from left to right, top to bottom, represent the example stars S0-2/S2, SO-102, Ea, Eb, 
Ec, and Ed, respectively. Open and solid circles mark the locations of the apoapsides and periapsides, respectively. Three total orbits for 
each example star are shown here and the first orbit of each star starts at t = 0. 


measurements that may be used to constrain the quadru¬ 
ple irroment and test GR. The cardioid-like shapes of the 
curves shown in Figure [8] are due to the position change 
of the star with respect to the central MBH, and con¬ 
sequently the change of the MBH potential aird the de¬ 
flection angle. If the true position of the star is right in 
between the MBH and the observer, the apparent posi¬ 
tion should be at the coordinate origin in Figure H] We 
note that the trajectories of apparent positions of a dis¬ 
tant star deflected by Jupiter, which is moving around 
sun, have similar car dioid-like shapes due to sim ilar un¬ 
derlying physics Isee iKooeikin fc Makarovir2007D . 

5.2. Spin-induced Redshift Differences 

Figure [9] shows the redshift evolution of each example 
star monitored by a distant observer for three full orbits. 
As shown in this figure, the maximum spin-induced red¬ 
shift (or blueshift) differences over a full orbit are about 
0.4, 0.05, 2.28, 82, 55, and 3.7kms-i for S0-2/S2, SO- 
102, Ea, Eb, Ec, and Ed, respectively. 

The redshift of a star is approximately given by 

Z ~ -V,. ~ -Vk,.' + |(Vk/c) 2 + crg/r + JV(a), (32) 

with the first, second, third, and fourth terms in the 
right-hand side representing the Newtonian Doppler shift 
due to the Keplerian motion of the star, the special rel¬ 
ativity correction to the Doppler shift (or the transverse 
Doppler shift), the gravitational redshift, and the part 
of redshift due to the MBH spin, respectively. Here the 


redshift Z is defined to have the same unit (kms ^) as 
that of velocity. In the Newtonian approximation. 


Vk,z' = sin/'[eorbCOsT'-bcos(T'-br;')], 

V «orb(l - e^^J 

(33) 

Vk (2GM./r - GM./aorb)^^" , (34) 

r ~ aorb(l - eorb)/(l -b Corb cos v'), (35) 

and (5V(a) is negligible comparing with those other terms. 
The second and third terms on the right side of Equa¬ 
tion (15^ are much less than the first term when r >> Cg. 
According to Equation (l32ll . approaches the maxi¬ 
mum redshift and the maximum blueshift (or the max¬ 
imum blueshift and the maximum redshift) at v' ~ T' 
and TT — T' if sin /' cos T' > 0 (or if sin /' cos T' < 0), re¬ 
spectively. At v' = T' (or v' = -K — T'), the star crosses 
the sky plane {x'y'). As seen from Figure [HI the approx¬ 
imations given by Equation (I32[) are well consistent with 
the full GR calculations. 

The gravitomagnetic held generated by the spinning 
MBH influences not only the motion of a star close to 
it but also the propagatio n of photons from the star to 
the distant observer (i.e.. lThornelll98CT . Both of these 
effects are encoded in the amount of shifts (or the red¬ 
shift) of lines in the star spectrum measured at different 
times0 Figure [TO] shows the redshift difference {5Z) of 

® The low-order effects, such as the GR gravitational redshift, the 
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Fig. 6.— Top panel: projected apoapsis displacements due to the 
MBH spin after the motion of one full orbit. The solid (dashed) 
line represents the results of apoapsis shift defined in the same way 
as ^apo for those stars having the same initial orbital elements as 
those of Ec (Ea), but with different eccentricities; the spin and its 
direction of the MBH is fixed at a = 0.99 and {i, e) = (45°, 0°). 
Dotted lines represent th e an alytical estimates of (5apo for those 
stars by using Equation inii. Bottom panel: solid line (dashed 
line) represents the difference of <5 between two apoapsides of two 
stars with the same initial orbital elements as those of Ec (Ea), with 
one of the two rotating around a rapidly spinning MBH (a = 0.99) 
and the other around a non-spinning MBH (a = 0). Here 5 is 
defined in the same way as that shown in Figure (O at apoapsis. 
The dotted lines are the same as those in the top panel. The 
star symbols mark the location of the star Ea, Eb, Ec, and Ed, 
respectively. 

each example star at each given moment for the case 
with a rapidly spinning MBH (a = 0.99) and that with a 
non-spinning MBH. For each star, as seen from this hg- 
ure, the redshift difference changes over the course of a 
single orbit, and it is most significant near the periapsis 
and most insignificant near the apoapsis; the amplitude 
of the differences also increases with increasing number 
of the orbital periods passed. For S0-2/S2, the redshift 
difference is about ^ 0.3km s“^ at the periapsis of its 
first orbit, and it mounts up to ~ 1km s“^ after about 
three full orbits; for the Eb and Ec, the redshift differ¬ 
ence can be up to > 100km s“^ at the periapsis of the 
third orbit. The amplitude of the redshift differences 
is the largest for Eb and the smallest for SO-102, which 
is fully consistent with the dependence of the redshift 
difference due to the Lense-Thirring precession and the 

relativistic Doppler shift, and the Rpmer delay, on the redshift have 
been investigated in t he literature separately through perturba tive 
approximations (e.g.. fucker et ^120061 : 1 Angelil Sahall2QlQI) . In 
this study, we only focus on the spin related line shifts (or redshift) 
and do not separately study those low-order effects one by one as 
in previous studies. The low-order effects, as well as the high-order 
effects, are automatically involved in our full GR calculations. 



t (yr) 


Fig. 7.— Dependence of the spin-induced position displacements 
on the spin direction and the star orbital orientation. The top panel 
shows the distance at any given moment between the apparent 
position of a star rotating around an MBH with a = 0.99 and 
that of the star around an MBH with a = 0, and these two stars 
have the same initial orbital elements as those of Ea. The solid, 
dashed, and dotted lines represent the cases for the MBH with a 
spin direction of {i, e) = (71°, 90°), (45°, 180°), and (159°, 63°), 
respectively. The bottom panel shows the distance at any given 
moment between the apparent position of a star rotating around 
an MBH with a = 0.99 and that of a star rotating around an MBH 
with a = 0; and these two stars have the same semimajor axis, 
eccentricity, and argument of periapsis as those of Ea, but with 
different inclination {!') and longitude of the ascending node (H'). 
The solid, dashed, and dotted lines represent for the cases with 
r) = (22°, 90°), (45°, 0°), and (73°, 62°), respectively. The 
spin direction of the spinning MBH is fixed at (z, e) = (45°, 0°). 
Three total orbits are shown here and the first orbit starts at f = 0. 

frame drag ging effect on the dista nce ( r) to the MBH 
revealed bv lKannan fc Sali^ (j2009li and lAngelil fc Sahal 
(|2010f) . i.e., SZ oc [and oc a~b(l ~ Sorb)”^ at the 
periapsis]. If the accuracy of the redshift measurements 
can reach < lkms“^ (or < 50kms“^), the spin of the 
MBH may be well constrained by using an observational 
redshift curve of S0-2/S2 (or the example star Eb or Ec) 
over three or more orbits and correspondingly > 45 years 
(or > 6 years for Eb, or > 1 years for Ec). 

The Lense-Thirring precession and the frame dragging 
can lead to a change in the orbital orientation of a star, 
as shown above and in Appendix [BJ One significant com¬ 
ponent of the redshift difference 6V{a) must be that in¬ 
troduced by the change of the orbital orientation, i.e.. 


(5Z(u)~ - . /- „ —r {cos/'[eorb cosT' -I- cos(T' -I- u')] 

y aorb(l ~ ^orb) 

xSI" + sin/'[—CorbSinT' — sin(T' -|- t>')]5T"} , (36) 
where 51” and 6T” are the changes of I' and T' induced 
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Fig. 8. — Differences in the apparent position of each example star at each given moment for the case with, and without, the consideration 
of the bending of light (obtained from the ray tracing technique). Panels from left to right, top to bottom, represent the example stars 
S0-2/S2, SO-102, Ea, Eb, Ec, and Ed, respectively. The open and solid circles mark the locations of the apoapsides and periapsides, 
respectively. Three total orbits for each example star are shown in this figure, but they are difficult to be distinguished. 


by the Lense-Thirring effect since the start of the orbital 
motion from Vq to v' (or Eq to E from tg to t). Further¬ 
more, 

SI" ~ [{E - Eg) - eorb(sinF; - sinF;o)]^/727r, (37) 

~ [{E - Eg) - eorb(sinF; - sinF;o)]^T727r, (38) 

Eg — CorbSinFio = —3—(^0 ~ Eg), (39) 

V “orb 


E-eoThSmE= J^-^{t-Tg). (40) 

V %rb 

Here SI' and JT' are the spin-induced changes of /' and 
T' per orbit (see Appendix B), Tg is the time of peri¬ 
centric passage with respect to the year of 2000, and 
tg = 2020 year is the starting time. The dashed lines 
in Figure [TUI represent the analytical estimates obtained 
from Equation (1361) , which seem to be roughly consistent 
with the numerical results (dashed lines) obtained from 
the full GR calculations, at least at the apoapsides and 
the periapsides. According to Equation (1551) and Ap¬ 
pendix (B] dZ is roughly proportional to a“j.| 3 (l — Corb)”^ 
[or a“j.^(l—Corb)”^] at periapsis [or apoapsis], if 1—Corb ^ 
1. According to Equation (I32p . the maximum redshift (or 
blueshift) difference induced by the MBH spin at periap¬ 
sis are 0.2 (-0.06), 0.04 (-0.02), 0.9 (-2.0), 28 (-68), 13 
(—28) and l.Okms”^ (—1.4kms“^) over a full orbit for 
S0-2/S2, S0-102, Ea, Eb, Ec, and Ed, respectively, which 
are roughly consistent with the full GR numerical results 


shown in Figure 1101 There are indeed some differences 
between the approximations given by Equation (155)) and 
those obtained from the full GR calculations, which are 
due to (1) the adoption of the mean change rates of fl', 
T', and I' (for the purpose of illustration here, we do 
not intend to adopt the evolution forms of those change 
rates as functions of v')] and (2) the neglect of the light 
propagation effects and high-order precession in the ap¬ 
proximations and the slight differences in the estimates 
of 5Z (as discussed for the position displacements above; 
see Fig. |5]). 

Figure ITT] shows the maximum redshift difference (top 
panel) and the maximum blueshift difference (bottom 
panel) induced by the MBH spin {SZ) obtained for stars 
with the same orbital elements as those of Ec (or Ea) ex¬ 
cept with various initial eccentricities. With an increas¬ 
ing number of orbits, the absolute value of the maximum 
redshift (or blueshift) difference mounts up. As seen from 
this figure, the analytical approximations are quite con¬ 
sistent with the numerical results, especially when the 
semimajor axis of the star is larger, which suggests that 
the main contribution to the redshift difference SZ is 
the change of the star orbital orientation and thus the 
change of the projection of the velocity to the line of 
sight. For stars with small semimajor axes, the analyti¬ 
cal estimates deviate from the full GR numerical results, 
which strengthens the necessity of using the full GR cal¬ 
culations for accurately constraining the MBH spin. 

Note that the redshift difference due to the spin ef¬ 
fects on the propagation of photons is about 10 ~ 100 
times smaller than that on the star orbit (see also 
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Fig. 9.— Redshift evolution curves of example stars. Panels from left to right, top to bottom, represent the example stars S0-2/S2, 
SO-102, Ea, Eb, Ec, and Ed, respectively. In each panel, the solid (open) circles mark the peric ent er (apocenter) passage of the star, and 
the solid and dashed lines represent the numerical results and the approximation given by eq. respectively. Totally three full orbits 

are shown here and the first orbit starts at t = 0. 


lAngelil fc Sahal 120101) . For S0-2/S2, if the accuracy of 
redshift measurements is on the order of > 1km s“^, 
the spin effect on photon propagation may be negligi¬ 
ble (although the Schwarzschild effect on photon prop¬ 
agation is sig nificant, close to 10km s~^, and cannot be 
neglected; see lAngelil fc Sahal 1201011 . However, for stars 
with smaller semimajor axes and higher eccentricities, 
such as the example stars Eb and Ec, it is necessary to 
include the spin effect on photon propagation (on the or¬ 
der of ~ 3km s“^ for Eb and ~ 1kms“^ for Ec, close to 
10% of the total spin-induced redshift differences) in or¬ 
der to get an accurate constraint on the MBH spin. Our 
full GR calculations are fast and efficient, which enable 
us to automatically have the high-order precession effects 
due the MBH spin and thus provide a more accurate fit¬ 
ting scheme for future high-precision measurements (see 
Section [6]) compared with those ap proaches adopting t he 
perturbative approximations le.g.. lAngelil et al]l2010t) . 

The spin-induced redshift difference {SZ) also depends 
on the spin direction and the orbital orientation of the 
star as shown in Eigure [T^l similar to that for the spin- 
induced position displacement shown in Eigure [H The 
top panel of FigurefT^ illustrates the dependence of 6Z on 
the spin direction by showing the evolution of dZ of Ea, 
for three different spin directions. The bottom panel of 
Figure [12] illustrates the dependence of 6Z on the orbital 
orientation of stars, by showing the evolution of 6Z of 
Ea and two other stars with the same semimajor axis 
and eccentricity as those of Ea but with different orbital 
orientations. As seen from this figure, these dependences 


help to simultaneously constrain both the absolute spin 
value and the spin direction. 

6. A GENERAL METHOD TO CONSTRAIN THE MBH SPIN 
PARAMETER 

In this section, we introduce a general method to con¬ 
strain the spin parameter and the metric of the GC MBH 
via the relativistic motion of a star in the immediate 
vicinity of the MBH by using the MCMC technique. 

We first generate mock observations for each example 
star listed in Table [1] in order to investigate whether 
the spin parameter of the GG MBH can be tightly con¬ 
strained by monitoring the orbital motion of a close-in 
star through the next generation facilities, such as the 
GRAVITY on VLTI, TMT, and E-ELT. For this pur¬ 
pose, we assume that the MBH mass is M, = 4 x 10^ Mq 
and the distance from the Sun to the GC is Rqc = 8 kpc. 
The MBH spin and its direction are assumed as those 
listed in Tabled) With each given set of initial param¬ 
eters for the six orbital elements of a star, we can ob¬ 
tain the trajectory of the star’s apparent position on the 
sky plane (ciobs.i, Pohsj) and the redshift curve {Zohsj) 
for any given time period through the full general rela¬ 
tivistic calculations described in Sections |3| and |4| Here 
the subscript j denotes a mock observation made at the 
observer’s time Gbs.i- The spin-induced position dis¬ 
placement is most significant near the apoapsis and the 
spin-induced redshift difference is most significant near 
the periapsis. The best strategy would be monitoring 
the star near its apoapsis and periapsis more frequently 
than at other locations in order to efficiently measure 
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Fig. 10.— Difference between the redshift of each example star at each given moment for the case with a rapidly spinning MBH {a = 0.99) 
and that with a non-spinning MBH (a = 0). Panels from left to right, top to bottom, represent the example star S0-2/S2, SO-102, Ea, Eb, 
Ec, and Ed, respectively. Open and solid circles mark the locations of the apoapsides and p eriap sides, respectively. Solid lines and dashed 
lines represent the numerical results and analytical approximations obtained from Equation I36[ respectively. Totally three orbits for each 
example star are shown here and the first orbit of each star starts at t = 0. 


the spin-induced effects and constrain the spin. Never¬ 
theless, here we set the intervals between two consec¬ 
utive observational times as dtobs oc ® (in the ob¬ 
server’s rest frame). For each star, the total number of 
mock observations is 120. Further assuming a set of ex¬ 
pected accuracies on the apparent position (cTp) and the 
redshift (cr^) measurements of a star, e.g., (cTp, az) = 
(10/ras, lkms“^), we obtain the mock observations on 
the apparent position and the redshift of the star. The 
points in the top two panels of Figures [13] and [T5|) show 
the mock observations of the apparent positions (top left) 
and the redshifts (top right) for S0-2/S2 and Eb, respec¬ 
tively. With those mock data, we can use the following 
procedures to constrain the MBH properties (mass and 
spin) and the distance from the MBH to sun. 

6.1. MCMC Fitting Procedures 

1. For a given set of parameters for the MBH 
spin parameters (a, i, e), the MBH mass (M,), 
and the initial orbital elements of a model star 
(aorb,o, eorb,o, ^ 0 , Iq, Tf), Up), the orbital motion 
of the star is generated over a period time that 
is the same as that of the mock one from full GR 
calculations as described in Sections El Note that 
the motion of the star is defined with respect to the 
central MBH, which is approximated to be at rest 
in the Galactocentric rest frame. 0 

® The acceleration of the MBH with respect to the observer are 
negligible on the timescale considered in this paper (< 50 years). 


2. According to the motion of the model star, we 
can obtain the apparent position (ctj, /3j) and the 
redshift Zo^sj of the star at each observation time 
tohs,j, for any given distance of Rgc from the MBH 
to the Sun (see Section |4|). In order to determine 
both the trajectory of the apparent positions of a 

First, the acceleration of the Sun with respect to the GC is on the 
order of 10~^knis“^ per year, and the total change of the relative 
velocity of the MBH due to this acceleration on the timescale < 50 

years is on the order of 10~^kms“^. Second, the Brownian motion 
of the MBH due to background stars is estimated to be on the or¬ 
der of 0.1k m s~^ and the chang e timescale may be larger than 10^ 
years (e.g.. [M^erritt et al.11200'^ . so that the velocity change due to 
the Brownian motion should be < 0.01kms“^ on a timescale of 50 
years. Therefore, the initial position of the MBH and the relative 
motion of the observer with respect to the MBH (to Sgr A*) can 
be approximately described by the following six parameters: the 
distance from the MBH (or Sgr A*) to the observer, the initial 
position of the MBH (or Sgr A*) on the sky plane (two parame¬ 
ters, and the proper motion and the radial motion of the observer 
with respect to the MBH (or Sgr A*; three parameters). For an 
accurate reconstruction of the orbital configuration of a star, it is 
necessary to also include all of those parameters in the MCMC fit¬ 
ting procedures described below. For the purpose of demonstration 
purpose in the present paper, however, we assume that the MBH 
position is fixed (or the same as that of the Sgr A*) on the sky 
plane and both the proper motion and the radial motion of the 
distant observer with respect to the MBH (or Sgr A*) are zero in 
the Galactocentric rest frame. Therefore, in the following MCMC 
fitting procedures, the MBH position is described by only one pa¬ 
rameter, i.e., the distance from the distant observer to the MBH 
(or Sgr A*). With future high-quality observations (e.g., on the 
position and the motion of Sgr A*), it is possible to also constrain 
all of those six parameters, simultaneously, in practice. 
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Fig. 11.— Maximum redshift difference (top panel) and the maxi¬ 
mum blueshift difference (bottom panel) induced by the MBH spin. 
The MBH spin and its direction are assumed to be a = 0.99 and 
(f, e) = (45°, 180°). Solid and dashed lines represent the numeri¬ 
cal results obtained in the same way as those shown in Fig. llOl for 
stars with the same orbital elements as those of Ec (solid lines) 
and Ea (dashed lines), respectively, except with various eccentric¬ 
ities. Dotted line s represent the estimates obtained approximately 
by Equation ilMt . This figure indicates that the spin-induced red- 
shift is mainly contributed by the change of the orbital orientation 
due to the Lense-Thirring precession and the frame dragging. 


star on the observer’s sky plane and the redshift 
curve, totally 11 total parameters, i.e., © = 

(-^GCi -^•5 G C, norb,0; (^orh.O: ^0’ ’ "^0’ ^o)? 

are needed. The exact values of the last six 
parameters for the initial orbital elements of the 
star are not interesting for the purpose of this 
study. 


3. To check whether the model star can fit the mock 
observations, we use the y^-statistics for both the 
trajectory of the apparent position (aj, f3j) and the 
redshift curve {Zj), i.e.. 


Xp 


E 


(Oj 0:obs,j) T (/3j /5obs,j)^ 


and 


xl 


= E 

i=i 


[Zj Zohs,jY 


(41) 


(42) 


respectively, where N is the total number of ob¬ 
servations. If both the trajectory of the apparent 
positions and the redshift curve of the mock star 
are used for the fitting, we have = Xp + Xzi 
however, if only one of the two data sets is used, 
we should have = Xp or x%- 


Fig. 12. — Dependence of the spin-induced redshift difference 
on the spin direction (top panel) and the star orbital orientation 
(bottom panel). Legend similar to Figure 0 

4. For a given set of mock observations D of the ap¬ 
parent positions, the redshifts, or both, of a star, 
the posterior probability distributions of those pa¬ 
rameters 0 can be constrained according to the 
Bayesian theorem, i.e., P(0|D) oc P(D|0)P(0) oc 
exp(—x^/2)P(0), where ^*(0) is the prior distri¬ 
bution of 0 and assumed to be flat over an ini¬ 
tial guess range for each parameter. The MCMC 
method is adopted, with an implement of the 
Metropolis-Hasting algorithm, in order to effi¬ 
ciently get the best fit to the data D and obtain 
constraints on the model parameters 0 by search¬ 
ing the parameter space. 

6.2. Fitting Results 

Figure [TS] illustrates the best fits to the mock obser¬ 
vations of S0-2/S2, including both the apparent posi¬ 
tions (solid line in the top-left panel) and redshifts (solid 
line in the top-right panel) The residuals of the fits are 
shown in the two middle panels. The solid circles in 
the two bottom panels show the difference between the 
mock observations and those predicted from a star with 
the same initial orbital elements around a non-spinning 
MBH, and the solid lines represent the differences be¬ 
tween the best fits and those predicted from the star 
around a non-spinning MBH, similar to 6 and SZ shown 
in Figures [S] and [TUI respectively. Obviously, the mock 
observations (points) are well fit by the model. 

FigurejTTjshows the two-dimensional contours and one¬ 
dimensional probability distributions of the parameters 
M,, i?GCj a I b and e, respectively. The parameters of 
the best fit are listed in Table [2] As seen from Figure (TTl 
and Table [H the absolute value of the MBH spin can 
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TABLE 2 

Best-fit parameters for the MBH mass and spin and the distance to the GC. 


Example Star 

a 

i 

n 

n 

<5Hgc“ 

(pc) 

(5M. » 
(lO®Af0) 

S0-2/S2 

0.812+°-i«/ 

O/ _27 

1567^1 

—0 c)Q+u.y» 
U.OU_Q 97 

-o.4i/]:;^ 

SO-102 

0-457tSifr 

09+78 

19+m 

-0.3l/];®® 

0.52/2;« 

Ea 

0.922tri27 

Cl +13 
^-*--13 

178/“ 

-0.34/];®® 

U.4D_2.o4 

Eb 

o 

00 

CD 

1 + 

o o 

b b 

(-• o 

ibtl 

ISl/j 

-1.16/y® 

-| t:Q+3.55 

j..OO_3 29 

Ec 

0.987tO;[](2 


178/g 

-1 07+^-®^ 
-*-•^'-2.61 

-0.80/®;" 

Ed 

0.840tO;ig 

6+23 

180/®® 

— 1 54+2®^ 

1 1 9+2.67 
-I--L^_2.79 


Note. — Notes. The error associated with each parameter represents the 2-a error of the best-fit value of the parameter. “ The best 
fit of the distance to the GC is Rqc = 8 kpc -|- (SRgc- ^ The best fit of the MBH mass is M. = 4 X 10® Mq -|- 6M,. 


be reasonably constrained by using the orbital motion of 
S0-2/S2 over two or three full orbits (30 — 45 years), if 
(1) the astrometric and the redshift accuracies can reach 
(Tp = 10/tas and az = 1kms“^; (2) the MBH spin a is 
close to one; and (3) the spin is pointing toward some di¬ 
rections with moderate or even maximum spin-induced 
effects on the trajectory of apparent position and the 
redshift curve, e.g., (i, e) = (45°, 200°), (49°, 126°) or 
(131°, 306°). According to Figures [5] and ITOl d mounts 
up with increasing number of orbits, therefore, the longer 
the time of monitoring the motion of a star, the tighter 
the constraints on the MBH spin obtained for a given 
set of (cTp, (Tz)- If the astrometric and redshift accu¬ 
racies can be a factor of a few times, or more, higher, 
i.e., dp < 3/ias and/or az < 0.3kms“^, the MBH spin 
can be constrained by using the motion of S0-2/S2 over 
only one full orbit (a period of ^ 15 year). However, 
if the astrometric and the redshift accuracies are a fac¬ 
tor of two times or more lower, i.e., dp > 20/ias and/or 
tr^ > 2kms“^, then the MBH spin cannot be constrained 
by using the orbital motion of S0-2/S2 within a few or¬ 
bits. 

Note here that both the distance to the GC and the 
MBH mass can be constrained simultaneously to an un¬ 
precedented high accuracy, i.e., 6M,/M, ^ a few times of 
10“^ and SRqc/Rgg ~ a few times of 10“"^. Even if the 
astrometric and redshift measurements are significantly 
less accurate, e.g., dp = 50/ias and/or dz = 10kms“^, 
the MBH mass and the distance to the GC can be still 
constrained to an accuracy of 0.98% and 0.36%, respec¬ 
tively, although the MBH spin cannot be constrained, by 
monitoring the motion of S0-2/S2 over less than three full 
or bits. These resul t s are consistent with those obtained 
bv lWeinberg et all (|2005H . 

If the absolute value of the MBH spin is smaller than 
0.99 by a factor of C, then it requires the astrometric 
and/or redshift accuracies of dp < lOC^^^as and/or 
tr^ < C^^kms”^ or a longer period for monitoring the 
star (> (2 — S)^ full orbits), in order to distinguish the 
MBH from a Schwarzschild MBH. 

According to Figures [7] and [HI the amplitudes of <5 
and SZ (due to the spin effects) also depend on the spin 
direction. However, the pattern of the dependence of 
S on (i, e) is somewhat different from that of SZ. For 
example, the amplitude of S is the largest (or moder¬ 
ate) when (i, e) = (72°, 90°) [or (45°, 180°)]; while the 
amplitude of SZ is moderate but not the largest (or mod- 



Fig. 13. — Mock trajectory of the apparent positions and the 
redshift curve of S0-2/S2 over three full orbital periods and their 
best fits. Top two panels are for the mock apparent positions (solid 
circles in the left panel) and the redshift curve (solid circles in 
the right panel) and their best fits (solid lines); middle panels are 
for the residuals of the best fits; and bottom panels represent the 
distances (solid circles in the left panel) or the redshift difference 
(solid circles in the right panel) between the mock observational 
positions (or redshifts) of S0-2/S2 and the positions (or redshifts) 
of an assumed star with the same initial orbital elements as the 
best fit of S0-2/S2 but around a non-spinning MBH, and the solid 
lines therein represent the differences between the predicted orbit 
of the best fit to S0-2/S2 and that of the assumed star. 

erate but close to the largest) when (i, e) = (72°, 90°) 
[or (159°, 63°)]. If both measurements on redshifts and 
apparent positions with compatible accuracies are avail¬ 
able, it would improve the constraint on the MBH spin 
and its direction. If only measurements on apparent posi¬ 
tions with accuracy of dp = 10/ras are available, then the 
spin is easier (or harder) to constrain by using the mo¬ 
tion of S0-2/S2 within several orbits if ( /, e) = (72°, 90° ) 
[or (159°, 63°)] (see more discussions in lYu et 5ll2015fl . 

Figure [12] illustrates the best fits to the mock ob¬ 
servations of Eb, similar to that shown in Figure [13] 
for S0-2/S2. The two-dimensional contours and one¬ 
dimensional probability distribution of the parameters 
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Fig. 14.— Two-dimensional probability contours and one¬ 
dimensional probability distributions of the best-fit parameters for 
S0-2/S2. For the two-dimensional probability panels, the color 
contour maps represent the mean likelihood of the MCMC sample, 
and the line contours represent the marginalized distribution, with 
the red solid line and the yellow dashed line as the l-cr and 2-a 
confidence level, respectively. The ‘-I-’ symbol in each color map 
represents the original set of the parameter value to produce the 
mock observations. The white dashed line in the panel for SRqq 
versus SM, represents M, oc (and 5M,/M, = ^SRqq/Rqq) 
with F ~ 2.4. For those panels showing the one-dimensional prob¬ 
ability distributions, the blue solid line and the green dashed line 
represent the one-dimension marginalized distribution and the one¬ 
dimensional mean likelihood, respectively. 
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Fig. 15.— Legend similar to that for Fig. I13l but for the example 
star Eb. 


Fig. 16.— Legends similar to those for Fig. 1141 but for the ex¬ 
ample star Eb. The dashed line in the panel for SRqc versus 5M, 
represents M, oc with F ~ 3.0. 

are shown in Figure 1161 and the best-fit parameters are 
listed in Tabled Since Eb has a smaller semimajor axis 
(300 AU) and a larger eccentricity (0.98) compared with 
S0-2/S2, the spin-induced effects on 6 and 6Z for Eb are 
larger than those for S0-2/S2 by a factor of several tens to 
a hundred, respectively (see Figs. ISl and ITOl). Therefore, 
the MBH spin and its direction can be more accurately 
constrained by using the motion of Eb over three full 
orbits (^ 10 years) if (ctp, tr^) = (lO^as, lkms“^), and 
the accuracy can be improved by more than one order of 
magnitude (see Tab. [2]). For much lower astrometric and 
redshift precisions, e.g., Up ^ 50^as and az ~ 10kms“^, 
which can be easily achieved by the next generation ex¬ 
tremely large telescopes, it is still possible to constrain 
the MBH spin with considerable accuracy by using the 
motion of a star like Eb over a period less than a few full 
orbits (< 10 years). 

To produce the mock observations, we have assumed 
that the MBH spin a = 0.99. If a is set to be smaller than 
0.99 by a factor of D, then the resulting amplitudes of S 
and SZ would be smaller than those shown in Figures [5] 
and [TUI by the same factor (see eas. lBSIlBSl andlUUI). If D is 
less than 10 and a > 0.1, we may still be able to constrain 
a or distinguish it from a Schwarzschild BH by using the 
motion of a star like Eb within a few orbits provided 
(cp, Tz) = (lO^as, lkms“^). If a < 0.5 and D > 2, 
however, it is difficult to get an accurate constraint on 
the spin value or distinguish it from a Schwarzschild BH 
by using the orbital motion of S0-2/S2 over a few or¬ 
bits, even if the spin vector is close to some specific di¬ 
rection [e.g., (45°, 200°) or (45°, 135°)] with which the 
spin-induced effects on 6 and/or SZ are most significant. 
In this case, a constraint on whether the GC MBH is an 
extremely rotating MBH or not is possible. If the astro- 
metric and redshift measurements are much more accu¬ 
rate (e.g., dp < 10D~^^as and/or crz < Il~^kms“^) as 
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discussed above, an accurate constraint on the MBH spin 
may still be obtained and it can be distinguished from a 
Schwarzschild BH. 

Figure [13 presents the one-dimensional probability dis¬ 
tributions of the best-fit parameters for other stars, i.e., 
SO-102, Ea, Ec, and Ed, respectively. The best-fit pa¬ 
rameters are listed in Table [3 As seen from Figure |T3 
and Table [U the MBH spin cannot be constrained by 
using the relativistic orbital motion of SO-102 within a 
period of less than two to three orbits (< 30 — 45 years) 
for tTp > 10 //as and uz > 1 km s“^, even if the spin di¬ 
rection is close to the one with the largest effects on 8 
and 8Z\ it can be well constrained by using other stars, 
i.e., Ea, Ec, or Ed. The reason is that SO-102 has the 
largest pericenter distance and thus the least changes of 
8 and 8Z per orbit caused by the spin effects. The MBH 
spin is less well constrained by using Ed, compared with 
those using Ea and Ec, because the pericenter distance 
of Ed is larger, compared to Ea, Eb, and Ec, though its 
semimajor axis is smaller than that of Ea and Eb. 

7. DISCUSSIONS 

7.1. Perturbations from Other Stars or Stellar 
Remnants 

In our full general relativistic calculations, the possi¬ 
ble perturbation on the orbit and motion of a star due 
to other background stars existing within the star or¬ 
bit is ignored. By assuming sever al different di s tribu - 
tions of stars close to the GC MBH. [Merritt et al.l (|201Clll 
performed post-Newtonian simulations and showed that 
the stellar perturbation may obscure the signals due to 
frame dragging and quadruple mome nt for stars beyond 
~ 0.5 —0.2mpc (~ 100 —400 AU); and lSadeghian fc Willi 
(IMlh also analytically investigated the perturbing ef¬ 
fects of background stars and found that it may still be 
possible to detect relativistic precession of stars within 
a few tenth milliparsecs. Therefore, it may not be easy 
to use S0-2/S2 (or even a star like Ea or Eb) to con¬ 
strain the GC MBH spin even if the accuracies of posi¬ 
tion and redshift measurements are sufficiently high (e.g., 
(Tp < 10/xas, (Jz < lkms“^). However, on the one hand, 
the exact number of stars or stellar remnants existing 
in the vicinity of the central MBH is not clear. It is 
not impossible that the stars or stellar remnants exist¬ 
ing within the S0-2/S2 (or Ea/Eb) orbit are very rare 
and their perturbations on the orbit and motion of SO- 
2/S2 (or Ea/Eb) are small and do not affect the spin 
signals much. Furthermore, the orbit-averaged torques 
from background stars are approximately constant in 
magnitude over year- to decade-long timescales, it is pos¬ 
sible to disentangle the stellar perturbatio n effects from 
those due to the MBH spin as discussed in iMerritt et al.l 
()20ir)tl . On the other hand, it is an important step to first 
check whether the spin-induced effects on the motion of 
a star, as close to the GC MBH as those example stars, 
can be detected by ignoring the perturbations from other 
possibly existing stars within its orbit, with future facil¬ 
ities with sufficient astrometric and redshift accuracies. 

Even if the perturbations from background stars on a 
target star are not negligible, the observational data col¬ 
lected around the pericenter, where the perturbations are 
substantially less significant than those near the apocen- 
ter, may still be useful to constrain the MBH spin pa¬ 


rameters (see lAngelil fc Sahal l2014f) . By using the red¬ 
shift data near the pericenter, strong constraints on spin 
parameters may be obtained since the spin-induced red¬ 
shift differences are the most significant near the pericen¬ 
ter. However, the constraints on spin parameters may 
be much weaker if only the position data are used be¬ 
cause the spin-induced position displacements near the 
pericenter are substantially smaller than those near the 
apocenter. As a check, we perform additional MCMC 
fittings for S2/S0-2 and Eb by only using the mock ob¬ 
servations near the pericenter (i.e., r < Oorb or r <C Oorb)- 
Our results show that the MBH spin parameters are 
poorly constrained by using the mock observations near 
the pericenter of S2/S0-2 (0 < a < 1, i = 61°t6io, 

and e = 160 °/j 2 j 3 o) because the spin-induced maximum 
redshift difference can only be marginally detected (see 
Figure [To]). For Eb, however, the spin parameters can 

be well constrained (a = 0 . 99 to'oi 2 i * = 45°tio, and 

I Q ^ X 

e = 182 °_ 3 o) because the spin-induced redshift differ¬ 
ences near the pericenter are substantially higher than 
the detection limit and can be accurately measured. 

It is also possible that an intermediate mass black hole 
(IMBH) exists in the v icinity of the GG MBH (e.g. , 
Yu fc Tremaind 120031 : iHansen fc Milosavlievid 120031 
Genzel et al. insir and this IMBH m ay perturb the or¬ 
bit an d motion of a star near the MBH. KIualandris et al.l 
(j 20 ir)l l have suggested that future observations on the 
S0-2/S2 orbit may be used to probe or constrain the ex¬ 
istence of an IMBH in the GC. 

We also note here that other effects, such as those from 
gas dynamics, wind mass loss, and tidal dissipation, may 
affect the orb it an d motion of a star c lose to the MBH. 
iPsaltisI (I20I2I1 and iPsaltis et al.l ()20I3r i have shown that 
these effects may be negligible and do not preclude the 
measurements of the MBH spin and quadruple moment. 


7.2. Accurately Constraining the MBH Mass and the 
Distance from the Sun to the GC 

With the expected astrometric and redshift accu¬ 
racies of the next generation facilities [(up, az) = 
(lO^as, lkms“^) or (50^as, 5kms“^)], the MBH mass 
and the distance from the Sun to the GC can be con¬ 
strained with an unprecedented accuracy on the order of 
a few 0.01% to 0.1 % or a few 0.1% to 1% if using S0-2/S2 
(see e arly works bv I Salim' fc Gouldlll999HWeinberg et al.l 
l2005f) . Even with such high accuracies of cTp and az, 
the degeneracy between the MBH mass and the dis¬ 
tance from the Sun to the GC remains, i.e., M, oc 
{8M,/M, = F,5i?Gc/i?Gc) and F ~ 2.4 - 3 0 (e.g.. see 
Figure s [U and flHll . similar to that found in iGhez et al.) 
(|2008l see Fig. 11 therein) and iGillessen et al.l (I2009L s^ 
Fig. 15 therein). Although the MBH spin can be better 
constrained by using the orbital motion of a star with 
a small semimajor axis and a high eccentricity, such as 
Eb, the constraints on M, and Rqc are not significantly 
improved. The tightness of the constraints on M, and 
Rgc depends on the ratio of the astrometric accuracy to 
the semimajor axis and the ratio of the redshift accuracy 
to the orbital velocity, but does not depend directly on 
the semimajor axis and eccentricity of the test star. 
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Fi g. 1 7.— One-dimensional probability distributions of the best-fit parameters for other stars obtained in the same way as done in Figs. ll4l 
and 1161 Top to bottom rows are for SO-102, Ea, Ec, and Ed, respectively. The solid line and the dashed line represent the one-dimension 
marginalized distribution and the one-dimensional mean likelihood, respectively. 


7.3. Usefulness of the Observations of Stars with Very 
High Eccentricities during Their Pericentric 
Passages 

In the GC, some stars may move around the MBH on 
nearly parabolic orbits (corb > 0.95) as suggested by the 
existence of the star S0-16/S14 ( with Cnrh = 0.963 an d 
pericenter distance of ~ 70 AU; iGillessen et al.l 1^00911 . 
With such high eccentricities, they may approach the 
MBH at a distance even substantially smaller than those 
stars with smaller semimajor axes, such as S0-2/S2, and 
thus more significantly affected by the spin effects. As 
demonstrated in Section[5l the spin-induced position dis¬ 
placements and redshift differences are most significant 
near the apoapsis and periapsis, respectively. Since the 
orbital periods of those stars on nearly parabolic orbits 
may be substantially longer than those with smaller semi¬ 
major axes (e.g., S0-2/S2 and SO-102), it is difficult to 
monitor them for a whole orbit. However, it should be 
useful to monitor the motion of those stars during their 
pericentric passage, when the spin-induced redshift dif¬ 
ferences (but not the position displacements) are most 
significant. We further check the validity of using the 
observations near the pericentric passage of a star on 
nearly parabolic orbit and find that the MBH may be 
accurately constrained if the pericenter distance of the 
star is smaller than a few tens AU. If there are more 
than one such stars, which would pass their periapsides, 
in the near future, monitoring their motions would be 
helpful in constraining the MBH spin. 

7.4. Direction of the MBH Spin 

The spin-induced GR effects, e.g., the Lense-Thirring 
precession and the frame dragging, on the orbit of a star 
depend on the direction of the spin. For a star with a 
given orbital orientation, the spin-induced effects on the 
changes of its apparent positions (or redshift) may vary 


by a factor of several to more than an order of mag¬ 
nitude if choosing significantly different spin directions, 
and the spin-induced effects are negligible for some spin 
directions, even if the absolute value of the spin is large 
and the semimajor axis of the star is small. However, 
the patterns of the changes of the position displacements 
and those of the redshift differences are different (see 
Figures [7] and [T^. Therefore, simultaneously using both 
the trajectory of the apparent positions and the redshift 
curve would help to obtain a better constraint on the 
spin, provided that the maximum values of the spin- 
induced position displacements and redshift differences 
over the observing period are substantially larger than 
the astromet ric and redshift accuracy, respectively (see 
discussions in lYu et ^1201511 . In this case, if only choos¬ 
ing one of the two data sets, the constraint on the spin 
would become losser because of the above orientation de¬ 
pendence. If the relativistic orbital motions of more than 
one star with significant different orbital orientations are 
available, the constraints on the MBH spin would also be 
tightened. 

7.5. Test for Alternative Metrics 

In this study, the metric of the massive object (pre¬ 
sumably an MBH) in the GG is assumed to be de¬ 
scribed by a Kerr metric. In principle, this object could 
be desc ribed by a metric, different from the Kerr met- 
ric (see [Willi 119 9311. e.g.. t he Johannsen-Psaltis metric 
l|Johannsen fc Psaltisll20lTI) . Similarly, it is also possible 
to use the relativistic orbits of stars in the immediate 
vicinity of the central object to check whether the cen¬ 
tral object is described by such alternative metrics or 
not. To do this, a modification to the numerical calcula¬ 
tions of the orbits and the ray tracing in such a metric 
is required, which is deferred to a future study. 

8. CONCLUSIONS 
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The innermost S-star discovered in the vicinity of the 
GC MBH is believed to be able to provide strong con¬ 
straints on the MBH spin and metric, and thus provide 
a strong test to the no-hair theorem of BHs and GR. In 
this paper, we develop a full general relativistic method 
to simultaneously constrain the MBH mass, spin, and its 
direction, and the distance from the Sun to the GG, by 
using the relativistic orbital motion of a star close to the 
MBH. 

The spin information is encoded in the trajectory of the 
apparent position of the star in the observer’s sky plane 
and the redshift curve measured by the observer, which 
include: (1) the changes in the orbital motion induced by 
the MBH spin and (2) the changes of the photon trajec¬ 
tory propagating from the star to a distant observer due 
to the spin. We numerically solve the motion equations in 
the Kerr metric for a star rotating around an MBH; and 
we follow the propagation of photons from the star to the 
distant observer by using a full general relativistic ray¬ 
tracing technique. Armed with these numerical methods, 
we can obtain the trajectory of the apparent position and 
the redshift curve of any star, with any given initial con¬ 
ditions, rotating around an MBH, with any given mass 
and spin (and the spin direction). We adopt six example 
stars with different semi-major axes and eccentricities. 
Two of these stars have orbital configurations similar to 
that of S0-2/S2 and SO-102, respectively; and the other 
four are adopted according to the probability distribution 
of the orbit of a star that is expected to be the closest one 
to the GC MBH. We obtain the trajectories of the ap¬ 
parent positions and redshift curves for those stars over a 
given period, and further obtain the mock observations 
of both the trajectories of their apparent position and 
their redshift curves by assuming astrometric (cTp) and 
redshift accuracies (ctz) of future facilities. According to 
our calculations, we emphasize some points below. 

First, the spin-induced effects on the orbital motion 
depend on both the absolute spin value (a) and spin 
direction. The amplitude of the spin-induced position 
displacements and redshift differences can differ signifi¬ 
cantly (by an order of magnitude) for MBHs with differ¬ 
ent spin directions but the same absolute spin value. 

Second, for S0-2/S2, the spin-induced effects would 
lead to a maximum change of 13/ras in the apparent 
position of the star within a full orbit for a maximum 
spin with direction of (49°, 125°) [or (131°, 305°)]. These 
value s are consistent with previous estimates ('e.g.-lwiil 
120081 : lAngelil et al.l 120101 1 . The maximum change de¬ 
pends on the spin direction. The changes in the apparent 
positions are mainly due to the spin-induced changes of 
the orbital motion, while the part due to the spin-induced 
changes of the photon propagation path is on the order 
of 0.01/ias, only ^ 0.1% of the total changes. 

Third, for S0-2/S2, the spin-induced effects would also 
lead to a maximum change of the redshifts (~ 0.3kms“^) 
over a full orbit for a maximum spin with direction of 
(129°, 171°) ]or (51°, 351°)]. The dependence of the max¬ 
imum changes in the redshift curves on the spin direction 
is different from that of the maximum changes in the ap¬ 
parent positions. The spin-induced maximum change of 
the redshift over a full orbit is also dominated by the 
orbital change, and the contribution from the photon 
propagation is on the order of 0.003kms“^, about the 


percentage level of the total changes. 

Fourth, for a star much closer to the MBH than SO- 
2/S2, such as Eb with a semimajor axis of ~ 300 AU 
and an eccentricity of > 0.98, the spin-induced effects 
would lead to a maximum change of 308^as in its ap¬ 
parent position within a full orbit for a maximum spin 
with a direction of ~ (111°, 275°) [or ~ (69°, 95°)]; and 
it would also lead to a maximum change of ^ 93km s“^ 
in the redshifts over a full orbit for a maximum spin with 
direction of ^ (150°, 319°) [or ^ (30°, 139°)]. The spin 
signals of a star like Eb can be more than an order of 
magnitude higher than those of S0-2/S2, and thus mon¬ 
itoring the motion of such a star could provide much 
tighter constraints on the MBH spin and metric. The 
maximum value of the spin-induced redshift differences 
over a full orbit is still dominated by the orbital change, 
and the contribution from the photon propagation is on 
the level of 10% of the total changes. 

Finally, in order to accurately constrain the spin and 
metric of the GG MBH (on the percentage level), it is 
necessary to consider both the spin effects on the stellar 
orbit and the photon propagation from the star to the 
distant observer. Although the changes in the apparent 
position of a star due to the spin effects on photon prop¬ 
agations are only about 0.001% — 0.1% of those caused 
by the change of the orbital motion, the changes in the 
redshift of a star due to the spin effect can be up to 
^ 1% — 10% of those caused by the change of the orbital 
motion if the semimajor axis is in the range from 1000 
to 100 AU. Therefore, it is important to include the full 
GR calculations of the photon propagations in order to 
accurately constrain the MBH spin. 

Adopting the MGMG method, we illustrate that the 
spin of the GG MBH may be able to be constrained by 
using the orbital motion of S0-2/S2 over a period of ~ 45 
years if (1) the spin value a is close to 1, (2) Cp ^ 10/ras 
and az ~ 1km s“^, which will be achieved by the new in¬ 
strument GRAVITY on VLTI, and (3) the spin is point¬ 
ing toward some specific direction, with which the spin- 
induced effects are likely moderate to most significant. 
If the astrometric and spectroscopic precisions can be 
further improved by a factor of several, then the GG 
MBH spin may be well constrained by using the rela¬ 
tivistic motion of S0-2/S2 monitored over a period of no 
more than 15 years. In the mean time, the distance from 
the Sun to the GG and the MBH mass can also be con¬ 
strained to an unprecedented accuracy (O.OI — 0.1%). If 
there exists a star with a semimajor axis of less than 
several hundred AU and an eccentricity of > 0.9 as ex¬ 
pected Isee iZhang et al.l 1201311 . the MBH spin can be 
constrained with high accuracy over a period of ~ 15 
years for cTp ^ 10/ras and az ^ lkms“^, even if a is only 
moderately large. The next generation facilities, such as 
the GRAVITY on VLTI, TMT, and E-ELT, may reach 
an astrometric accuracy of 10 — 50/ras[] Those extremely 
large telescopes may reveal some stars, if existing, with 
semimajor axes less than a few hundred AU and eccen¬ 
tricities > 0.9, within the S0-2/S2 (or SO-102) orbit. The 
spin (and its direction) of the GG MBH is likely to be 
tightly constrained by using the relativistic motions of 
such stars. Long time monitoring of the orbital motions 

’’’ http://www.tmt.org; http://www.eso.org/sci/facilities/eelt/ 
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of stars in the vicinity of the GC MBH by the next gener¬ 
ation facilities is likely to provide an ultimate dynamical 
test for the first time to the no-hair theorem and GR. 
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APPENDIX 


A. PHOTON ORBITS IN A KERR SPACETIME 

The motion of a photon in a Kerr spacetime is given by Equations ([3|)-([6|). We assume that an observer, located 
at {ro,9o,(j)o) = (10®rg,z,0) in the BH rest frame, receives photons emitted from a star located at R* = (r*, 0*, 0*), 
which is rotating around the BH. Define u = cos 6, the motions of thos e photons in the r — 9 plane are given by 
(|Bardeen et al.l[l97l : iCunningham fc Bardeenl[l973l: iChandrasekharlfTQSSU 



dfj, 


(Al) 
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where fj-o = cosi, = cos0*, and rhit —^ r* is the redshift of a photon when ^ = /i*, i.e., the photon hits the surface 
of the star at In Equation (lAll) . the left side is integrated from a distant observer to the surface of the star 

and it should be negative as tq > Chit ■ We should carefully choose the sign in front of the integration at the right side, 
i.e., if /io < /i*, t he si gn should be and if /tq > /i*, the sign should be “+”. Below, we denote the left and right 
side of Equation (EH) as and ^r, respectively. 


A.l. Integral over the 0-direction 

We denote the integral over /i on the left side of Equation m as 

dji 




vW)’ 


where 


Q[P) = q^ + (a^ -q^- - (1^), 

1 


= 


2a2 L 


y/ (A^ + =F (A^ + - a^) 


(A2) 

(A3) 

(A4) 


and —n+ and /r+ are the two real roots of 0(/i) = 0. Because 0(/i) must be > 0, the physical allowed region for a 
photon is < /I < M+' Once fi reaches the turning point (either /i_|_ or —^+) and bouncing back, the /i-integral at 
the right side of Equation (IA5|) changes the sign. If the integration range of the 0-integral is from //(> 0) to we 
may denote the function ^^(/x) explicitly as 




rv+ 


dp, 


If —/x+ < /i < 0, we have 


dp 


rt^+ 


dp 


M \/0(m) 

f^+ dp 




(A5) 


(A6) 


•Am \/0(m) Jo y/QiP) J-v \/0(Ai) 

According to iBvrd fc Eriednim] (jl954 p48, eq. 213.00), =^^(/x) can be expressed in terms of an inverse Jocabian 

= 5/xCn"^(cos(/5|m^) = ((/3|m^) 


elliptic function 
where 


9^^ = 


1^91 + 9 - 


rrin = 




M — 9 , 9 ’ 

Mi + M- 


(A7) 

(AS) 

(A9) 


cn“^(cos ip\mfj,) is the inverse Jacobian elliptic function, and E((/3|m^) is the incomplete elliptic integral of the first kind. 
Therefore, the integral ^/^(m) obtained from the integration over any range of ^ can be expressed as a combination of 
inverse Jocabian elliptic functions. 

For the purpose of tracing photons from a distant observer back to the surface of a star located at (r,^, 0*, (/>*), we 
integrate the trajectory of each photon from /xq to ^*. For those photons that can finally reach the star surface, some 
of them may reach the star surface when they first approach /x = /x*. For these photons, the /x-integral can be 
obtained by 

-[^^(/xo) - ^^(/x*)], if/3 > 0 and /xq < /X*, 

-[=^/x(Mo) + =^/x(m*)]: if /3 > 0 and /xq > /x*, 

+ [=^/x(/xo) - =^/x(m*)]) if/3 < 0 and /xo > m*, 

-k[-2^^(0) -k =^^(Aio) + if /? < 0 and /xq < /x*- 




^,1 


(AIO) 


Some photons may not reach the star surface when they first reach /x = /x*. However, they may reach the star surface 
after their trajectory, turning back from /x_|_ or —/x+ to /x* again. For this case, we have 



^IJ.,2 = 


,even5 
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Consequently, some photons may reach the star surface after their trajectories turn back to /i* for the third time. For 
this case, we have 



Li,2 

(A13) 

where 

f -2[2^^(0) - if /3 > 0 and /xq < 



1 -2[^^(0) - ^^(/x*)], if/3 > 0 and /xq > AX*, 

1 - 2 [^^( 0 ) - ^^(/x*)], if/3 < 0 and AXo > A^*, 

[ -2[2^^(0) - ^^(ax*)], if /3 < 0 and axq < A^*- 

(AM) 

Therefore, if some photons reach the star surface when their trajectories reach u* for the 2fc—th time [fc(> 1) is a 

integer], then we have 

u,2k — ^ 1 “1“ event 

(A15) 

while we have 

^ Li,2fc+1 — ^ Li,2fc ^^^,oddt 

(A16) 

if some photons reach the star surface when their trajectories reach /x* for the (2k + 1)- 

-th time. For most photons 


included in this paper, fc < 2 — 3. 


We denote the integral over r as 


A.2. Integral Over the r-direction 


^r{r) 



dr 

\/W)’ 


which can also be expressed in the form of the inverse Jacobian elliptic functions. To do so, we need to first solve 
equation R{r) = 0 and find its four roots, denoted as ri, r 2 , ra, and r^, respectively. Since R{r —>■ ±oo) —>■ oo and 
R{r = 0) = the solution of R{r) = 0 can be divided into two cases: case A, Rjr) = 0 has two real roots and 

two complex roots; and case B, R{r) = 0 has four real roots fe.g..see rChandrasekhai][T9^ iR.auch fc Blandfordlll9^ 
ICadez et ahlfTMl lYu fc L1[20^: iL^lT^iroi: lird~m[2nM IDexter fc AgollMIl . Below, we consider these two 
cases separately. 

Case A: we assume that ri and r 2 represent the two complex roots, and r^ and r 4 (r 4 < r^) represent the two real 
roots, respectively. Therefore, we have ri = where represents the complex conjugate of r 2 . Since 'rir 2 r 3 r 4 = 
= —a^q^ < 0, r 4 must be < 0 and r^ > 0, but r^ cannot be equal to 7 - 4 . Since R{r) must be > 0, we have 
either r > r^ or r < r^, and then the physically allowed region for r must be r > ra. We assume that the real part of 
the two complex roots ri and ra is u, while the imaginary parts of ri and r 2 are w and —re, respectively. Therefore, 
the two complex roots can be expressed as ri = u + iw and r 2 = u — iw, respectively. The two real roots can be 
expressed as 7-3 = —u + v and r 4 = —u — v, respectively, where u is a real positive number. Now we can denote the 
function ^rijr) explicitly as 


■^r2{r) = 


dr 


dr 




y/{f - ri){f - r 2 ){f - r 3 ){r - r 4 ) 


According to IBvrd &: FriedmanI (|1954 pl35, eq. 260.00), the above function can be expressed in terms of an inverse 
Jacobian elliptic function, i.e., 

where 


(fi = cos 


-1 


= 52Cn ^(cOS(A>|7772) = 52A((A>|7772), 

(A17) 

Vab’ 

(A18) 

(A + 5)2 - 47;2 

^2 = -:- 5 

4AB 

(A19) 

A = ^/{v — 2?x)2 + w'^, 

(A20) 

B = y/ (v + 2?x)2 + w'^ 

(A21) 

(A — B)r + (?x + v)A — (u — v)B 

(A22) 

_ (A + B)r + (?x + v)A + (tx — v)B 

cos^ ip <1. Therefore, the elliptica' 

integral F{ip\m 2 ) in Equation (IA17I) 


is well defined for r > r^. Therefore, may also be expressed as a combination of inverse Jacobian elliptic functions. 
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In this case, some photons may reach the star surface with r monotonically decreasing from tq to rhit, and some 
other photons may reach the star surface with r first decreasing from rg to rg and then turning back to rhit- Hence, 
we should have 

= -[^r2{'ro) - ^r2(?'hit)], (A23) 

and 

+ ,^r-2(?'hit)], (A24) 

respectively. 

Case B: assuming that ri > r 2 > rg > r^. Since rir 2 r 3 r 4 = —a^q^ < 0 and ri + r 2 + ra + r 4 = 0, we have at 
least r 4 < 0 and ri > 0. Therefore, the allowed physical regions for photons with R{r) > 0 are r > ri (region I) and 
rs < r < r 2 (region II, if r^, >0). For those photons in region I, i.e., r > ri, we denote the function ^rif) as 


■^r4.i{r) = 


df 


VW) 


df 


(r - ri)(f - r 2 )(r - r 3 )(r - r 4 ) 


(A25) 


According to iBvrd fc FriedmanI (jl954 pl28, eq. 258.00), Equation (IA25I1 can be expressed in terms of an inverse 
Jacobian elliptic function, i.e., 

=^r4,i(r) = (;4sn“^(sin (^|m4) = g4F{ip\m4), 

where 


94 = 


Lp = sm 


\/(?’l - ?'3)(?'2 - Ti)' 

( I {r 2 - r4){r - ri) \ 

I V (^1 “ ’'4)(r - r2) / ’ 


(A26) 

(A27) 


(A28) 


(4-2 - r3)(ri - r4) 

(ri - r3)(r2 - r4) 


(A29) 


For this region, if ri ^ r 2 , 0 < 7714 < 1 and thus Equation (IA26I) is well defined. If ri = r 2 , TO 4 = 1 and sin^ ip = 1, and 
thus Equation (jA26p oc sn“^(l|l) = 00 , which is not well defined. If ri = r 2 , however, the integration in Equation (jA25l) 
can be directly obtained as 



df 

VW) 


’^r4,leq{r) 


1 

\/{ri - r^){ri - r^) 


In 


\/{r - r-i){r - rq) 

r — n 


rf + + 2rir 

(r - ri).^(ri - r 3 )(ri - r4) 


(A30) 


For photons in the region I, some photons may reach the star surface with r monotonically decreasing from rg to 
rjiit, while some other photons may reach the star surface with r first monotonically decreasing from rg to ri and then 
turning back to the star surface. If ri ^ r 2 , hence, we should also have 

= -[’^r4,l{ro) - ^r4,lirhit)], (A31) 

and 

J^r = -[=^r4,l(7’o) + ^7-4,1 (7"hit)], (A32) 

respectively; if ri = r 2 otherwise, we should have 

= -[^4,Ieq(?’o) - ^4,Ieq(?'hit)], (A33) 

and 

= -[^4,Ieq(cc) + J^^4.Ieq(?'hit)], (A34) 

respectively. 

For those photons in the region II, i.e., r 2 > r > rg, we denote as 


=^r4,Il(7’) 




df 

ri){f - r2)(r - r3)(f - r 4 )' 


(A35) 


According to IBvrd &: FriedmanI (|I954 pll6, eq. 255.00), Equation (IA35I) can be expressed in terms of an inverse 
Jacobian elliptic function, i.e., 

■^r4.uir) = g4sn"^(sin(^|m4) = g4F{p\m4), (A36) 


where 


(f = sin 


/ j in - r3){r2 - r) \ 
\\j [r 2 - r 3 )(ri - r) j ' 


(A37) 
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54 , and 1714 are the same as that given by Equations (IA27I) and (IA29I) . respectively. Similar to Equation (IA26|) . 
Equation (IA36I) is well defined if ri r 2 as 0 < m 4 < 1 and sin^ < 1. If ri = r 2 , then Equation (IA36I) can be 
directly replaced by 


dr 






\/(?'l - »'3)(»’l - ’’ 4 ) 


In 


\/(r-r3)(r-r4) 
Ti — r 


rf + r-3r4 + 2rir 


(ri - r)y/{ri - r3)(ri - r4) 


(A38) 


For photons in the region II, some of them may reach the star surface with r monotonically decreasing from rp to 
rhit, while some others may reach the star surface with r first monotonically decreasing from rp to rp and then turning 
back to the star surface. For these two cases, if rp < r 2 and r 2 r^, we should have 


r 4 ,11(^0) ^ r 4 ,Il(^hit)], 


and 


= —2^r4,Il(?’3) + [=^r4,Il(»’o) + =^r4.II(^hit)], 
respectively; if r-p < ri and ri = r 2 otherwise, we should have 

— ['^^4,IIeq(^o) ^^^^4,IIeq(^hit)], 

and 

= —2^4,Ii(?' 3) + [^4.IIeq(?'o) + ^4,IIeq(?’hit)], 


(A39) 

(A40) 

(A4I) 

(A42) 


respectively. 

For the problem considered in this paper, the photons that we are interested in must be transported from the stars 
rotating around the central MBH to a distant observer. Our calculations show that their trajectories should not be 
bended too much to oscillate within the region II (i.e., < r < r 2 ) for many times. 


A.3. Thit 

According to Equation ini, SO rhit can be derived from this equatio n once the two integrals are obtained 

as described above. If a photon reaches 5 = /i* for k times, then leas. IA15l and lAl^ . 

Case A of the r—integral: 


1. If a photon reaches the star surface before it encounters the turning point rp, then 



^ Li,A: — r2('Co) ^ r2 ('Chit,A:)] • 

(A43) 

Therefore, we have the solution 





= ■^r2^ id^fJ.,k + ^r2{ro)). 

(A44) 

We define S21 = ^fj.,k + ^r2('Cp), then we may express rhit,fc as 



61-62COSV? 

A’hit,fe — ^r2 (^) “ 5 

a2 COS Lp — ai 

(A45) 

where 





oi = A — B, 

(A46) 


(X2 — A-\- 

(A47) 


bi = -{u + v)A + ( 

u — v)B, 

(A48) 


62 = -(u + v)A - ( 

u — v)B, 

(A49) 

and 

f S21 

cos (/? = cn — 

V 92 

7712 ^ , 

(A50) 


and 1712 is given by Equation (jA19l) . 


2 . 


If a photon reaches the star surface after it encounters the turning point rp, according to Equation (jA24l) . then 
we have 


= ■^r2{s22) 


hi — 62 cos Lp 
a 2 cos (/? — ai ’ 


(A5I) 


where 

S 22 = —d^iL,k — ■^r 2 {ro), (A52) 

and oi, 02 , 5i, & 2 , and cos(/? are given by Equations (IA46I) - (IA50I) . If S 22 > =^r 2 (oo), then there is no solution for 
rhit.k, and therefore it is not necessary to go to larger k to search for the solution of rhit,fe. 
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Case B of the r-integral: 

1. For photons in the allowed region r > ri (but ri ^ r 2 )., if they reach the star surface before they encounter the 


turning point ri, according to Equation (I A3II) we should have 

_ f _r2C-ri 

— ^r-4,l(®4l) — (J _ ^ ' 

'S 4 I — ^ [i,fc “1“ ^r4,l(r'o); 

- f4 . 2 

C = -sin If, 


where 


4-2 - 4-4 

S 41 

sin (^ = sn - 

54 


7714 


(A53) 

(A54) 

(A55) 

(A56) 


(74 and 7714 are the same as that given by Equations (IA27|) and (jA29l) . respectively. 


2. Eor photons in the allowed region r > ri (but ri ^ 72 ), if they reach the star surface after they encounter the 
turning point ri, according to Equation (IA32I) we should have 


— '^^ 44 ( 542 ) — 


where 


720 - 7l 

C-l ^ 

^42 — u,k ^ r4,l(4’o)5 


, S 42 

Sin (^ = sn - 

54 


7714 


(A57) 

(A58) 

(A59) 


C, 54 , and 1714 are the same as that given by Equations (IA56|) . (IA27I) and (IA29I) . respectively. If S 42 > —^^,k — 
^r 4 .i{oo), then there is no solution for rhit.fc, and it is not necessary to go to larger k to search for the solution 
of 7\,i\^ f„. 

3. Eor photons in the allowed region r > ri and ri = 72 , if the reach the star surface before they encounter the 
turning point ri, according to Equation (IA33I) we should have 


where 


4’hit,fc — '^^4,Ieq('®43)) 
543 ~ ^u,k “t” '^^4,Ieq(4'o)- 


(A60) 

(A61) 


4. Eor photons in the allowed region 7 > 74 and ri = 72 , if they reach the star surface after they encounter the 
turning point 74 , according to Equation (IA33I) we should have 


where 


7hit,k — '^^^4,Ieq('^44)5 
^44 — ^u,k '^^^4,Ieq(4'o)- 


(A62) 

(A63) 


If S44 > —J^fj,,k — =^4,11(00), then there is no solution for rhit,fe, and consequently it is not necessary to go to 
larger k to search for the solution of rhit,fe. 

5. Eor photons in the allowed region 72 > 7 > 73 , they cannot transport to the distant observer if 7q ^ 72 , and 
therefore are not interested in this study. 

A.4. Longitude and Time 

The motion of a photon in the (j)- and t-directions can be also obtained as (e.g., see lRauch fc Blandfordl[T994f) 


Ar^ + 2r(a — A) dr 


‘ Sign 


and 


t = r. 


- 2 r + 

^ + 2ar{a — A) dr 


sign 


— 2r o? 


VW) 


^sign 


+ 0. 


sign 


dfi 

(A64) 

2 2 

(A65) 
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provi ded th at th e mot ion of the photon in the r- and 0-directions are given (see Sections lA.ll and IA.2I) . In Equa¬ 
tions and (IA65|) . rgign and 0sign are the s igns of the r- and 0-integrals, which can be obtain by so lving Equa- 

tion (IA1|) (see details in Sections IaT 2I and IA.1|) . The integrations on th e ri ght sides of Equations (IA64|) and (IA65|) 
over either r or 0 can be expanded similar to those done in Sections IA.2I and lA.ll according to rgi gn an d 0sig n and how 
many times that they change sign; we do not show it here in detail. The r-terms in Equation s TA64D and (IA65I) can 
be accurately integrated by adopting the Gauss-Kronrod integration scheme (see also lBauch fc Blandfordlll994l l. The 
0 -terms can be obtained analytically by the elliptical functions as 






dpi 


rA‘+ 


dp 




dp 




-1 M 

--cos — 

1 - M+ V d'+ 


n^, rrifj, ) - (cos ^ ^ 


according to lBvrd fc Friednranl (11954 p48, eq. 213.02), and 




Jf^(/r)= / P 


dp 


rt^+ 


|a| 4 \ d+- m' 


:dfl - 




M ’ 


dfl 


= {g+ + (^cos 


-1 _ 0 _ 
d'+ 


- fj-_g^3^ cos 


A* \/0A*(m) 
-1 d 


d+ 


m. 


(A66) 


(A67) 


according to IBvrd fc Friednranl (|1954 p48, eq. 213.01), where S’, and 11 ar e the elliptical integral of the first, 
second, and third kind, respectively, g^ and are given by Equations (IA8I) and (IA9I) . and is given by 


= 



(A68) 


A.5. Apparent Position of a Star Close to the GC MBH in the Observer’s Sky Plane 

It is necessary to develop a fast and accurate method to find those rays propagating from the observer to the star. 
The most simple and direct method is to perform ray tracing calculations for a large number of rays with various (a, 
jS) and search for the minimum of the distance function d{a,/S) = (i(rhit, 0*, (/>hit) (defined by Eg. E71) . However, this 
method may be not efficient in finding the ray that actually hits on the star at any given position because (1) not every 
ray with arbitrary (a, /3) can end up at the surface of the cone with 0 = 0*; and (2) the contours for d{a, ft) appear to 
be extremely narrow elliptically in some cases; thus, the convergence for numerically finding the local minimum of d 
is slow. 

In this study, we implement a new method to find the local minimum of d{a,ft). We assume a flat space with 
coordinate {r',0',(f') as an approximation to the Kerr space. In this flat space, those rays from the distant observer 
that can hit the position (r',0*,(()') have 

a = r-'sin0* sini))', /3 = r'cos 0* sin0o — r'sin0* cos <()'cos0o. (A69) 

Here and hereafter the superscript ' indicates that a quantity is defined in the flat space. Since the example stars are 
still far away from the event horizon, the trajectories of photons emitted from those stars do not deviate away from 
the corresponding cases in a flat space too much. Therefore, the local minimum of the distance function d can be 
found by the following procedures. 


1. We first assume (rg, t/fg) = (r*, </>*) to obtain (ag, /3o), which is taken as the guess position of the star [(r*, 0*, 0*) 
in the Kerr space] on the observer’s sky plane. Adopting {a, ft) = {ao,Po)y we perform ray tracing calculations 
in the Kerr metric to hnd the position of the ray hitting the surface 0 = 0*, i.e., (rg, 0*, (()g), which may deviate 
from the real position of the star (r*, 0*, <()*) because of the curved spacetime. We calculate the distance between 
(r*,0*,(^*) and (rg,0*,(()g) as dg = d{ao,/3o). 

2. We select a number of points (e.g., eight points) on a circle with radius of dg = dg around (r-Q,0*,^Q). For 

the j-th point (j = 1,...,8) selected on the circle, i.e., (r( ^), we can obtain (aij,fiij) according to 

Equation (|A69p . Using {aij, Pi j), we perform ray tracing calculations again and find the position {rij, 9i,, pij) 
that the ray hits on the surface 9 = 9i,, and consequently, we obtain the distance of that position from the star 
as di,j = d{aij, Pij). Among those (eight) points on the circle, we find the one [denoted by (ri^m, d*, (^ii^m)] 
with the smallest distance d (denoted by di^m), and the corresponding point in the flat space is denoted by 

3. If dl.ni < C^* (f = 10“® ~ 10“®), then we find the apparent position of the star and the calculation is terminated; 

if c^i.m > and di^m > then we adopt the same point (rg,0*,(()g) in the flat space as the center of a new 
circle with a radius of d'^ = dg/5; otherwise, we choose the point (r[ 0*, (f[ jjj) in the flat space as the center of 

the new circle and its radius is set to be d[ = dg. 
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4. We select a number of points {r2j,&*,4’2,j) circle and obtain {a 2 ,j, I32,j)- We further perform ray 

tracing calculations and find the position {r2j,0^,(j)2j) that the ray hits on the surface 9 = 0*. We calculate 
the distance of that position from the star i.e., c? 2 j = d(a 2 j, P 2 ,j) in the Kerr space. Among those 

selected points, we find the one with the smallest d, i.e., d 2 ,m, and the corresponding point in the flat space is 
denoted as {r'2^rn^9^,4>2^^). 

5. Similar to the above procedures 3 and 4, we do iterations until we find the apparent position of the star at 
{ri,,9^,(f>^) in the observer’s sky plane. 

The above method is fast in searching for the apparent position (a, fd) of a star at (r*, 0*, (/>*) in the Kerr spacetime 
on the distant observer’s sky plane, which enables an efficient MCMC fitting to the mock data obtained for those 
example stars. 

B. ANALYTICAL ESTIMATION OF THE SHIFT OF STAR APPARENT POSITIONS DUE TO THE LENSE-THIRRING 

PRECESSION AND THE FRAME DRAGGING 

The Lense-Thirring effect causes the precession of the orbital plane via the change of the longitude of the ascending 
node at a mean rate of 

A 2 aM. 

and the change of the argument of periapsis at a mean rate 

< t >=-SficosC, (B2) 

where O and T are the longitude of the ascending node and the argument of periapsis defined for the stellar orbital plane 
with the equatorial plane of the BH as the reference plane, and C is the angle between the orbital angular momentum 
of the star and the spin direction and given by C = arccos(cos /' cosi +sin/' sini sin(f2' — e)). The accumulated changes 
of and T per orbit are 


3/2 


Sn = 27r , (B3) 


a, 


and 


3/2 


(5T = -27r- 


6 are 


3/2 




cosC. (B4) 

-‘orbV^ ~ 

Note that most of the changes in 12 and T take place near periapsis. The values of the changes are Sid = 0.0017°, 
0.00059°, 0.010°, 0.14°, 0.074°, 0.0092°, and (5T = 0.0044°, 0.0010°, -0.015°, -0.21°, -0.11°, and -0.014° for SO- 
2/S2, SO-102, Ea, Eb, Ec, and Ed, respectively. Considering the projection effect, the changes of it', T', and /' can be 
obtained by 

cos /' sin i sin(f2' — e) 


<512' = 


cosi — 


sin /' 


^ 12 , 


<5T' = - Si"^sm(12'-e) ^^ 
sin /' 


(5T, 


and 


61' = sinicos(12' — e)<512. 


(B5) 

(B6) 

(B7) 

—>■ 0 , 


If ^ 0°, i.e., the star orbit is close to being on the MBH equatorial plane, /'—?►/, 12' — e —?► 7 r/2 , then <512' 

(5T' —?► —2(512, and 61' —>■ 0. 

After the motion of one full orbit, the difference between the position of an example star in the sky plane at 
apoapsis/periapsis for the case with a rapidly spinning BH (a) and that with a non-spinning BH (a = 0) is 


(5apo/peri — aorb(l ± Corb) [i5^12'(l — sin^ T' sin^ /') -|- (5^T'(1 — cos^ T' sin^ /') -I- sin T'^ sin^ I'6^I' 
+2 cos I'Sit'ST' - 2 sin T' cos T' sin /' cos I'ST'ST - 2 sin T' cos T' sin TSil'ST] , 


(B8) 


where the signs ‘-I-’ and are for the difference at apoapsis ((5apo) and periapsis ((5peri), respectively. 

We may als o in clude the high-order precession due to the quadruple moment of the MBH by replacing <512 and <5T 
in Equations (IB3D and (IB4I) by <512 -|- <512q and (5T -|- (5Tq, respectively, where 


(512q = — cos^- 


37ra^rg 


"^orb 


( 1-0 


and 


5To = 


1-5 cos^ C Z'Ko^r'i 


orb 


( 1 - 0 ^' 


(B9) 


(BIO) 











See detailed derivations from [Barker fc O’ComieUl (|1975f) and iWex fc Kooeikinl (|1999ll . 










